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Measured Properties of Strong ‘ Unipotential ”’ 
Electron Lenses 
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MS. received 27th May 1948, and in amended form 18th October 1948; read at Science Meeting 
3rd December 1948 


ABSTRACT. ‘The Hartmann test and the “ knife-edge ’’ method for the determination 
of lens constants are described as applied to the investigation of electrostatic “‘ unipotential ’” 
lenses. ‘The results of measurements on a number of such lenses are presented graphically. 

It is shown that the curvature of the principal surfaces plays an important part and that it 
is, therefore, necessary to distinguish between longitudinal spherical aberration and difference 
of zonal focal lengths. Up to about half the lens diameter, primary spherical aberration 
only is noticeable in most cases.. The dimensionless spherical aberration constant S can be 
approximately represented, over a wide range of modifications of lens design or lens poten- 
tials, by the formula S= K,(f/D)?, f being the focal length and D the internal diameter of the 
centre electrode of the lens. K, is a constant of the order of 10 for the lenses investigated. 
without stops limiting the axial length of the field, and of the order of 20 for lenses incorpor- 
ating such stops. In the case of one specially constructed lens, good agreement is found with 
calculations by Ramberg. In an experimental unipotential lens from which the diverging 
parts of the electrostatic field were removed, the spherical aberration is as low as in magnetic 
lenses of comparable dimensions. 


sy INTRODUCTION 

HE “‘unipotential”’ or “‘single”’ lens occupies a unique position amongst 

the known electron lenses. It shares with the other electrostatic lenses the 

advantage that its properties are independent of the ratio e/m of charge to 
mass of the focused particle and depend only on the ratios of the lens dimensions 
and on the ratios of the voltages applied to its electrodes (provided that the particle 
speed is so low that no relativistic corrections are required). Its distinguishing 
feature is that its electrodes need only be connected to the cathode potential and to 
the final anode potential, which is mostly earthed. This simplifies the electrical 
circuits associated with the lens and makes the lens action largely independent of 
variations of the potential difference between cathode and anode. 

The two basic forms of the “‘ unipotential”’ lens are three perforated diaphragms 
with aligned circular holes, and a short cylinder placed between two longer 
cylinders. The outside electrodes are connected to the same potential, V,, and 
the centre electrode to the cathode of the electron-optical system, or to any other 
suitable potential V.AV,. If V,A0, the advantages of making the lens action 
independent of the lens voltages is lost, but greater flexibility in operation is gained. 
If V,, is sufficiently negative, the lens is converted into an electron mirror. 
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The “‘unipotential”’ lens is used in cathode-ray tubes, in electrostatic electron 
microscopes, and in ion focusing devices. The image-forming electron bundles 
are here frequently of such angular aperture that the aberrations of the focusing 
lens may become of considerable importance. Spherical aberration in particular 
limits the performance of cathode-ray tubes (Maloff and Epstein 1938, p. 132) and 
of electron microscopes (Rebsch 1938). Whilst the theory of the “‘ unipotential” 
lens has been discussed extensively (Scherzer 1933, 1936, Rebsch and Schneider 
1937, Plass 1942, Ramberg 1942, Bachman and Ramo 1943, Bruck and Romani 
1944, and others), little attention was given to the measurement of their properties 
except by Johannson and Scherzer (1933), Gobrecht (1941), and Mahl and 

Recknagel (1944). 
Because of this lack of information, the writer studied experimentally several 
such lenses during 1942-1944, mainly to find out whether modifications of the 
basic design could give low spherical aberration values. 


§2. GEOMETRIC-OPTICAL RELATIONS 

For the “‘ unipotential” lens, the refractive index of the medium is the same in 
image space and object space, and the lens formula can be written in the simple 
form 

aPrib=lif” eee (1) 
where a is the distance of the object P and b the distance of the image Q from their 
conjugate principal surfaces, and f is the focal length—see Figure 1, where the 
principal surfaces H and H’ are shown crossed-over and curved, as is found in 
actual measurements. All unipotential lenses investigated were of the symmetri- 
cal type; in this case the principal surfaces are symmetrical about the mid-plane M 
of the lens and intersect the optical axis at a distance d)/2 from the mid-plane, dy 
being the separation of the apexes of the principal surfaces. ‘The dimensions a, 
b and d, are not directly measurable; the values determined by the experiment are 
the object distance p and the image distance g from an arbitrary reference plane, _ 
e.g. mid-plane M, and the angles « and £ of the rays in the object and image spaces. 

If we let several rays from the object point P fall on the lens at different heights 
of incidence 7, (measured at R, the point of intersection of ray and principal 
surface H), we find that the image points Q,, are displaced in an axial direction from 
the paraxial image point Q, towards the lens. The difference gy —q,, = Aq, is the 
longitudinal spherical aberration. ‘The th ray cuts the plane perpendicular to the 
optical axis through the image point Qp, the ‘“‘ Gaussian”’ image plane, at a distance 
6, from the axis. 6, is the lateral spherical aberration (radius of the ‘disc of 
confusion”’ due to the mth lens zone). 

It is customary to refer the aberrations to the focal plane, i.e. to consider the 
aberrations produced by a bundle of rays entering the lens parallel to the optical 
axis. ‘The distance of the focus of the mth zone from the mid-plane will be 
denoted by g*. From the theory of lens aberrations it is known that the primary 
longitudinal and primary lateral spherical aberrations of the nth zone can be 
written in this form (fy = paraxial focal length) : 


BG Solty fo OS Se IE eee (2) 
and oF aa Sfol7n/fo)?- sees (3) 


S is a dimensionless constant characterizing the spherical aberration of the lens 
(Bachman and Ramo 1943), dependent on the lens geometry and relative refracting 
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power, but independent of the lens dimensions. It is therefore particularly 
convenient in a comparison of the aberrations of various types of lenses; it is 
related to the frequently used aberration constant Cx (Cosslett 1945) by the formula 


OBIS hi ee en (4) 


It is seen from equation (1) that for a finite object distance, as indicated in 
Figure 1, formulae (2) and (3) are replaced by Ag, =(bo/fo)(bn/f,)Ag* and 


Figure 1. Diagram of optical paths through (thick) lens. 


5,,=(r,,/b,)Ag,. Neglecting second-order terms, we obtain from these equations 
two alternative expressions for the constant S : 


De Sa ea) lo! Open Wi eRe a (5) 
and S= (8,,/o)(fol Tn)? Die Cae (6) 


An assumption, made by most authors, is that the principal surfaces are planes. 
‘The present measurements showed however that they are more nearly paraboloids. { 
Hence, a meridional section through each principal surface can be expressed in the 
form 4d,,=4d,—1r2(T —4)/fo. The shift Ad/2 of the principal surface from its 
mid-plane distance d)/2 on the axis is therefore given by the equation 


Ad,,/2 = aA 1 2) /fo- 5 Aaordto (7) 


The dimensionless constant T is proportional to the coma of the lens. A 
value T =0 would represent a slight curvature of the principal surface such that, in 
the absence of spherical aberration, S =0, the lens would fulfilthe ‘‘sine-condition” 
i.e. would give uniform magnification for off-axis points. 


+ Ramberg denotes the aberration constant S' by C, Scherzer and others Cg by Co. 


A curvature of the principal surfaces had already been observed in magnetic lenses (Becker and 
Wallraff 1938), but has not previously been found in electrostatic lenses. 
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If the focal length is measured, as usual in light optics, along the ray, we 
have for. the zonal difference of focal length Af, =f, —fo (provided po>/fo) : 
Af,, =Aq* —4Ad,, +12/2fo. Hence 

Mf, =iSs=lyiy a aa a ee (8) 
or S=T+7Af,,/7.. eee (9) 
This zonal difference of focal length produces a zonal difference of magnification 


AM,,, as M,,~f,,/Po for Po>fo. AM,, is therefore proportional to (S—T). 


$3, TEST iE PHO Ds 


Two test methods were employed, the ‘‘ Hartmann Test” and the “ Knife 
Edge Test”. 


(i) In the “Hartmann Test”? (Hartmann 1904), a well-known optical test 
method which has been applied by previous authors (Maloff and Epstein 1938, 
p. 130, Klemperer and Wright 1939) to electron lenses, the lens zones are deter- 
mined by placing a diaphragm with a number of pin-holes in front of the lens. 
The pencils thus selected are observed in the image space on a screen which can 
be moved to different positions (Hartmann’s method of extra-focal observation). 
The pencils can be located in the image space (angles f,, and intersection distances 
g,) to within 1%, whilst the position of the pencils in the object space (angles «,, 
and distance p, =p ) can be calculated with an accuracy closer than 0-3°% from 
the geometrical arrangement of electron source and zone selection disc. Dimen- 
sions are chosen to make py > q,,, in order to approach the case of parallel incidence. 

The lens constants f, and d,,, or fy, S and T, can then either be worked out from 
the zonal magnification M,, =tana,/tanf, and f,=M,p,)/(1+™,), which may 
be improved by corrections, or by a graphical method. In this, the incident and 
emergent rays are plotted on graph paper on a large scale, object point P and 
image points Q,, being placed on the same side of the lens. ‘The incident ray is 
plotted above and the emergent ray below the optical axis. Figure 2 represents. 
such a plot (lens No. 4) The intersections of the conjugate rays then define the- 
points R,, on the principal surface. Hence, d,,/2, q, and 6, can be taken from the 
graph. Then a, and b, are easily found, and f,, is obtained from the lens formula 
equation (1). 

To find S and T, the fractions Aq,,/r2, 5,,/73, 4d,,/r2, and Af,,/r2 are formed ; 
their average values are then inseried into equations (5), (6), (7) and (9). Thus, 
each value of S and T is obtained from a number of measured points, giving a 
fairly reliable result. Agreement of the three values of S, determined from rela- 
tions (5), (6) and (9), was usually within 10°% or better. ‘The S values given below 
are the averages of these three evaluations. 

The accuracy of the final values of fp is estimated as closer than +1%, of S 
about +10%, and of d)/2 and T about +20%. 


(ii) In the “ Knife Edge”? method the aberration is determined from the 
distortion of a straight edge, e.g. the edge of a narrow slit (see also Boersch 1939, 
Dosse 1941). It is particularly applicable to lenses which show little aberration, 
or when it is desirable to observe only for small angular apertures « and B. It 
was, however, not recognized by the earlier authors that the test as described by 
them did not give S but (S—T) and was, therefore, not a true test of spherical 
aberration ; only in bad lenses, as we shall see later, is S so much greater than T” 
that (S—T) can be considered as nearly equal to S. 
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If we place a “knife edge” a distance y from the axis (Figure 3) and illuminate it 
with a bundle of rays originating in P, we should, in the absence of spherical 
aberration, obtain a “shadowgraph”’ of the knife edge at the position B, on the 
screen, located at a distance / from the mid-plane M of the lens. 

As the ray after leaving the lens does not pass through the paraxial image 
point Q, but through the zonal image point Q, corresponding to the zone distance 
7, from the axis (7, being the off-axis distance of point R,, the projection of the knife 
edge on to the principal surface H from the point P), the “‘shadowgraph”’ will be 
displaced to the point B,. 


-40 


Figure 3. 


“200-80 -60.-40-20—0 
mm 


Figure 2. Experimental plot of Figure 3. Diagram of “ knite-edge = test: . 
electronpencils through lens Figure 4. Determination of displacement Xo of apex 
under test. of parabola in “‘knife-edge’’ test. 


From the geometrical relations of Figure 3 it can be derived (see Appendix I) 
that the displacement Ah=h,—hy of the point Q, from point Q, is given, for 
Afifo<1, by Ah =lr,Af/f5. 


In view of equation (8), this can be written as 
AUT Siscins swe aes (10) 


The shadow image of the straight ‘“‘ knife edge”’ will therefore appear curved, 
showing as a rule pincushion distortion. It can be proved from the cubic Jaw of 
primary spherical aberration that the distorted image of the edge is a parabola, the 
origin of which is displaced by the amount x)=/(S — T)(71/fo)’, the same value as 
given by equation (10). If the “knife edge” forms the edge of a slit placed 
symmetrically about the axis and a square is inscribed into the centre part of the 
“‘shadowgraph” image of the slit (see Figure 4), then the displacement of the 
image measured along the sides of the square equals 2x). Hence the required 
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displacement x, of the apex of the parabola can be easily obtained as the difference 
of the half-width of the square and the half-width of the “waist” of the shadow 
image; (.S—7) is then found from equation (10). f) can be determined from the 
angular magnification as discussed before, and / and r, are given by the geometry — 
of the test arrangement. 

The accuracy of measuring (8 — 7’) by this method can be greatly increased if 
the radius of curvature p of the parabola is used as a measure of the distortion. 
This can be done by placing stencils (e.g. circles drawn on a sheet of celluloid) on a 
photograph of the distorted shadow image of the “knife edge”. This method of 
measuring (.S — 7) is very accurate because in matching a stencil to a smooth curve 
the average value of a great number of measuring points is obtained. Values of 
Af/fo as small as 0:01 have thus been measured with an accuracy closer than 
LO. 

Under the condition p> fp, which can usually be arranged, (.S — 7) is then given 
by equation (11), which can be derived from the value of p as function of the 
parameters of the parabola : 


GSD )( bf) fol 2 aps Sf | eee (11) 

It should be noted that the ‘“‘ knife edge’’ method can be made to give S and T 
separately as well as fy and d,/2, if the distorted ‘“‘shadowgraph”’ is recorded for 
two different screen positions. If the distance of the knife edge from the axis is 
kept constant, these two screen positions have to be near the paraxial image point 
Q), preferably within a distance equal to about half the focal length on either side of 
Qo, to yield four independent equations for the four unknowns fo, dy/2, S and T, two 
equations being given by the magnifications and two by the distortions in the two 
screen positions. If both screen positions are on the same side of the focus and a 
fairly large distance from it, the four equations are no longer sufficiently indepen- 
dent to give the required accuracy in evaluating the four lens constants separately. 

This difficulty can be overcome by observing the magnifications and distor- 
tions of the knife edge for two (or more) different, sufficiently distant screen ~ 
positions for several off-axis distances of the knife edge (or for several slit widths), 
and then using the graphical method described in connection with Figure 2, in 
conjunction with the (S—T7) values obtained from the distortion measurements... 

These conditions were only fully recognized after the experimental work 
described in this paper was completed. As a constant slit width was used, only 
Jo; 4/2, and (S—T) were measured. 


84 EXPERIMENTAL TECHNIQUE 


Most lenses investigated were built up from equal diameter tubular elements, 
and were constructed as parts of small sealed-off cathode-ray tubes. 

The construction of such a “‘lens-test tube”’ is shown diagrammatically in 
Figure 5. The electron beam was produced by a triode electron gun of conven- 
tional type. The anodes A, and Aj, forming the outer elements of the unipoten- 
tial lens L, and its centre electrode A, were made from 4in. diameter seamless 
non-magnetic nickel copper alloy tubing (bore D=11-7mm.). The electrodes. 
were mounted in the usual way on mica discs held rigidly in position by four stout 
non-magnetic alloy rods. These rods extended about 4in. beyond the end of the 
electrode A, and served as guide rods for the sliding fluorescent screen S. The 
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whole assembly was mounted on a stem F and sealed into a tubular internally 
“metallized” envelope E, the end part of which consisted of a light walled bulb free 
from optical imperfections. 

A “zone disc”’ Z (or a parallel edged slit when the knife-edge method was used) 
was arranged, as a rule, in a field-free position within the anode cylinder A. Six 
pairs of holes at different off-axis distances were drilled in the zone disc, allowing 
the simultaneous observation of six zones. The positions of the holes of each 
zone disc were carefully measured with a projecting microscope. 

The fluorescent screen could be moved by tapping the tube. ‘The screen 
position could be read on a scale fixed to the outside of the tube envelope. The 
image of the electron pencils on the fluorescent screen was photographed in a 1:1 
ratio, and all measurements were taken on the photographic plate. The positions 
of the rays were measured for four screen positions, moving the screen from one 
extreme end tothe other. ‘To check the correctness of the scale readings of the 
screen positions, photographs were taken as well of the ‘“‘shadow images”’ of the 
zone disc or slit. "These shadow images were obtained as the straight line projec- 
tions from the point P of the holes in the zone disc when V,, was made equal to V,, 
i.e. if the lens action was removed (Figure 2). 
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Figure 5. Experimental “ lens test ’’ cathode-ray tube. 


The influence of the field of the four supporting rods which penetrates through 
the gaps between the centre electrode of the lens and its outer electrodes was found 
negligible in all but the largest lens zone used, becoming just detectable in the 
outermost zone (r,=D/4). 


$5 Nisa b Bis, UNV ose AT ED 

Fourteen different unipotential lenses were investigated. ‘The essential 
design parameters measured in lens diameters D, except for lenses Nos. 13 and 14, 
where the dimensions are given in millimetres, are shown diagrammatically in 
Figure 6. 

The lenses Nos. 1 to 6 are of the same type and may be arranged in ascending 
order of lens thickness, although a lengthening of the centre electrode, leaving the 
gaps between electrodes unchanged, is more effective in thickening the lens than 
leaving the centre electrode unchanged and lengthening the gaps. ‘The relative 
order of lenses Nos. 4 and 5 might therefore be reversed. 

Lens No. 7 is similar in dimensions to lens No. 3, but contains two diaphragms 
across the outer electrodes of the lens, each provided with acentre hole. ‘This lens 
was made to study the influence of the holes in the diaphragms on the spherical 
aberration. As it appeared that the field near these holes was still rather high, 
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lens No. 8 was constructed; this is similar to lens No. 7, but has the diaphragms 
removed further into the outer electrodes, thus placing them into a rather weaker 
field. 

Lens No. 9 was made in order to study the properties of an electrostatic lens 
from which the diverging field sections have been removed. It is identical in 
dimensions with lens No. 7, but the diaphragms are replaced by the zone disc on 
one side and by an open wire mesh disc (120 meshes per inch, about 70% trans- 
parency) on the other side. Although each hole in the zone disc acts as a lens, the 
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Figure 6. ‘Types of investigated lenses. 


principal ray passing through the centre of the hole is not deflected from its position, 
and the effect of such “‘ pinhole lens”’ is only a decrease or increase in the size of the 
cross-section of the electron bundles which pass through it and enter the main 
lens.* A similar argument can be applied to the section of the wiremesh. The 
field shape in this lens is such that it is purely converging, whereas all other 
unipotential lenses have diverging fields at the entrance to, and exit from, the lens. 


* This point has since been discussed more tully by Gabor in connection with ion-tocusing grids 
(Gabor 1947). ‘ 
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Hence lens No. 9 has rather interesting properties, although it is not a practical 
‘proposition as a focusing element in ordinary applications. Further, this lens 
allows a clear-cut comparison with lens No. 7, being only different from this by the 
-absence of the entrance and exit holes. 

Lens No. 10 also is almost identical with lens No. 7, but the ends of the outer 
lens electrodes which face the centre electrode are hemispherically curved, in an 
-attempt to change, at least in part, the curvature of the equipotential field surfaces 
near the outer lens electrodes in a more favourable direction. 

Lenses No. 11 and 12 as well are similar to lens No. 7, but incorporate conical 
inserts (cone angle 120°) which serve to reduce the field where the entrance and 
exit apertures are situated. The only difference between these two lenses is the 
‘presence of apertures in lens No. 12. 

Lens No. 13 is similar to lens No. 8, except for a rounding-off of the edges of the 
-electrodes. It was used as objective lens in an experimental three-stage electro- 
_ Static electron microscope constructed by the writer some years ago. This lens 

‘is not the best one tested from the point of view of spherical aberration, but was 
chosen because it appeared to lend itself to comparatively easy manufacture and 
-easy alignment, and because it showed a quick decay of the electric field outside the 
lens. 

Lens No. 14 corresponds to the lens ‘‘’T'ype B” described and analysed in its 
action by Ramberg (1942), and was made in order to see how closely theory and 
experiment would agree. 

Lenses Nos. 13 and 14 were tested by the knife-edge method, al] others by the 
Hartmann method. 

Sor LES Th RESULTS 

The measured lens constants as functions of relative zone radius 7/D are plotted 
in Figures 7 and 8 for lenses Nos. 1 to 12, all for the condition V,/V,=0. These 
figures give the relative focal length //D ,the relative focal distance g*/D from the 
mid-plane of the lens, and the relative distance $d/D of the principal surface from 
‘the mid-plane for rays of parallel incidence. It can be seen how necessary it is to 
distinguish between longitudinal spherical aberration Ag, =(q% —¢*) and zonal 
-difference in focal length Af, =(fo—/,)- 

In Figure 7 we see that a progressive increase in the lens thickness leads to an 
increase in refracting power; the longitudinal spherical aberration decreases in 
step with the paraxial focal length fp, so that Af/f) or Ag/gy remain almost constant, 
‘for a given off-axis distance 7/D. The parameter T remains practically constant 
within measuring accuracy, at a value around 6. 

This picture changes when we come to the very thick lenses, Nos. 5 and 6. 
‘The paraxial focal length of these lenses has not decreased further, and the curv- 
_ature of the principal surfaces has increased instead, becoming very great for lens 
No. 6,the parameter TJ rising to 10 in lens No. 5 and to 14 in lens No. 6. The rise 
‘in T is, indeed, so pronounced that the spherical aberration constant S, which 
showed a decrease with increasing lens strength in the first lenses, rises sharply, 
-although the term (S— 7) shows a slight further drop. Lens No. 6 is particularly 
‘interesting; the zonal focal length, which decreases for small off-axis distances, 
increases again for greater distances, becoming equal to the paraxial focal length fo 
at a distance 7,/D =0-2, even increasing beyond the paraxial value for greater values 
-ofr/D. 'Thus the curve of f/D resembles the familiar aberration curves of spheri- 
-cally corrected optical systems. Nevertheless, lens No. 6 shows a high degree of 
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spherical aberration, given by the g*/D curve, as the principal surfaces recede 
much more quickly from the mid-plane of the lens for increasing off-axis distances. 
than the zonal focal length increases; the shape of the f/D curve is mainly due to 
the fact that the paraxial focus lies within the lens field. This behaviour of electron 
lenses for large off-axis distances, which was found to a lesser degree also in some: 
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Figure 7. Relative focal length f/D, relative focal 
distance (for rays incident parallel to optical 
axis) g*/D, and relative separation of principal 
surfaces from mid-plane of lens $d/D as function 
of relative off-axis distance r/D for lenses 
Nos. 1-6. 
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8. Relative focal length f/D, relative focal’ 
distance (for rays incident parallel to 
optical axis) g*/D, and relative separation of 
principal surfaces from mid-plane of lens. 
$d/D as function of relative off-axis 
distance r/D for lenses Nos. 7-12: 


Figure 9.  Paraxial refractive power D/f, of 


lenses Nos. 1-6 as function of approximate: 
length of lens L/D. 


other cases (see for instance lens No. 1, Figure 7, and lens No. 12, Figure 8), may 
suggest a possibility of obtaining zonal correction. 
the writer that a lens such as No. 6 might be used with advantage as a partly 
distortion-free projector lens in electron microscopes. 


Mr. 'T. Mulvey pointed out to 
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In Figure 9 the relative refracting power D/f, for lenses Nos. 1-6 is plotted as. 
a function of lens thickness L. Whilst the lens thickness would be more correctly 
defined as the length over which the electrostatic field exceeds a certain minimum: 
value, it was taken here, for practical reasons, as the distance between the planes. 
bisecting the gaps between inner and outer lens electrodes. This definition, which 
does not take full account of the actual length of the refracting field, is probably 
the main reason for the scatter of the points in Figure 9. The refracting power 
D/fy appears at first to go up in proportion to the lens thickness L, as one would 
expect from the formulae for the focal length of a weak electron lens, but seems to 
approach a maximum value of the order 1-2 to 1-5 (for V,=0). 

Turning now to lens No. 7, we see that the effect of the perforated diaphragms. 
is a slight weakening of the refracting power (by 19°), whereas the spherical 
aberration constant S shows an increase by a factor of 2:8. The influence of the 
entrance and exit holes on the performance of the lens is thus very marked, contrary 
to the expectations of earlier authors (Bachman and Ramo 1943). This deterior- 
ation of performance is not altogether surprising owing to the addition theorem of 
the spherical aberrations of several lenses (see Appendix II), and the fact that the 
relative off-axis distance used by the electron beam passing through the entrance 
and exit holes is very much greater than in the main lens. 

Lens No. 8 shows a rather small relative focal length, indicating that the removal 
of the entrance and exit holes into regions of a weak field has been effective in 
eliminating the influence of these holes on the focal length; yet the increase of the 
spherical aberration over that of lenses without diaphragms is still appreciable. 

The properties of lens No. 10 fall half way between those of lenses Nos. 7 and 8,, 
as might be expected from its geometry. 

A comparison of lenses Nos. 11 and 12, which are similar to lens No. 8 in 
dimensions and optical constants, brings out once more very clearly the influence: 
of apertures not situated in an absolutely field-free region. ‘The paraxial focal 
length remains unaffected but the spherical aberration is doubled. 

We find a reverse effect in lens No.9. The elimination of the diverging parts 
of the lens field leads to a decrease in focal length and to a considerable improve- 
ment in spherical aberration, putting this lens into a class of its own amongst the 
electrostatic lenses investigated. The spherical aberration constant of this lens is 
very nearly the same as that of a magnetic lens of the same //D ratio. 

A summary of the measured lens constants is given inthe Table. ‘To bring out 
further the connection between focal length and spherical aberration and coma, the 
values of the aberration constants S and (S — 7) are plotted in Figures 10 and 11 as 
functions of relative focal length f)/D. Figure 10 indicates that for lenses whose 
basic design elements are cylinders the parameter S can be written as * 


Hee Re siT ese MS Gar) AV (12) 


The remarkable point is that the constant K, does not appear to depend, 
within wide limits, on the dimensions or even the shape of the lens elements, yet 
it is quite sensitive to the introduction of entrance and exit stops, as seen from the 
difference of the two curves in Figure 10. The very thick lenses Nos. 5 and 6 are 
omitted since they do not fall on the curve, the constant S being rather higher. In 
contrast to this, lens No. 9, the purely converging lens, has an aberration constant: 


* Lenses whose elements are diaphragms do not seem to follow such simple law. 
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which is smaller by a factor of about 2-5 than would be expected from the lower 
curve in Figure 10. The constant K, has a value of 9-5 for the lower curve, and a 
value of 18 for the upper curve. 

In Figure 11, the parameter (S— 7) is plotted for all lenses except No. 14. 
Apart from the very thick lenses, for which (S—T) falls away more quickly, the 
parameter (.S— 7) can be represented by the formula 


(S=D= KD 4. eee (13) 


Again the constant K, is about twice as great for lenses with entrance and exit 
stops as for the others, being about 10 for the former and about 5 for the latter. 


Table 
Lens fo/D $d,)/D 1S db S—T 
No: Ve/¥a (41%) (420%) (410%) (20%) (425%) arias 
1 0 3-10 0-11 91 6-0 85 
2 0 1:92 0:10 35 6-3 29 
3 0) 1539 0-09 16:5 4-8 IIE?) 
+ 0) 0:92 0-09 9-0 6-0 3-0 
5 0) 0-97 0-09 125 10-0 25) 
6 0) 0-86 0-11 it7-0 14-0 3-0 
7 0 1-65 0-04 46 7:4 39 Diaphragms 
8 0 0:87 0-04 14-2 6°61 29°017"6 Diaphragms 
9 0 0-78 0-09 2:6 2-0 0:6... Converging only 
10 0 INSU 0-06 23-5 10-0: SiS Diaphragms 
Mit 0 0:93 0-09 U5) 3°8 37 Conical inserts 
12 0 0-93 0-09 12-0 3-6 8-4 As No. 11 but 
(410%) diaphragms 
13 +0-20 2-43 (0-11) — = FI | 
+0-10 170 (0-11) — — 30 Similar to No. 8 
0 12? 0-11 — — 10:5 > but rounded 
—0:10 = 0°83 (0-11) — — 4-7 electrodes 
—0-20 3:58 (0-11) — —— Ds Badal 
14 +0-30 4°38 0-12 — — WE ) 
+-0-20 2:82 0-12 oo — 295 
+0:10 2-02 0-12 ss = 15-5 r Ramberg’s “ B ” lens 
0 1:38 0-12 — — 8-4 
—0-10 0-94 0-12 + —- Sy? 


J 


The full circles (intermediate curve) refer to lens No. 13, for which only (S — T) 
could be measured, the different values of f)/D being produced by changing the 
ratio V/V, of the electrode voltages. From this we conclude that, within wide 
limits, the aberration parameter (S — 7) depends mainly on the relative focal length 
Jo/D, whether this be brought about by a change in design or by an adjustment of 
electrode potentials, the values of S or (S — T) in different types of the investigated 
lenses not varying by more than a factor of two except in extreme cases. 

A comparison of theory and experiment are given by Figure 12, where the 
refracting power D/f, and the parameter (S— 7)D?/f? are plotted for lens No. 14 
(Ramberg’s “ B” lens) as a function of {(V,—V.)/V,}*. It has to be remembered 
that Ramberg’s calculated ¢/M value is not (S — T)D?/f2 but SD?/f?. Nevertheless, 
the agreement between Ramberg’s calculations and these measurements is very 
satisfactory, the constant 7 probably being rather small inthis case. This lens is 
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indeed better for a given f)/D value than the other lenses measured (except lens 
No. 9), but the difference is small, and still within a factor of two. 


Figure 10. Spherical aberration constant 
S as function of relative paraxial focal 
length f,/D. x, lenses Nos. 1-4; 
O, lenses Nos. 7, 8 and 10.- 
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Figure 11. Aberration term (S—T) as 


function of relative paraxial focal 
length f)/D for lenses Nos. 1-13. 


e Lens N° 14. 
ALens N°13. , 


Asymptote 


Wee 


Refractive power D/f and aberration constant (S—T)D?/f? as 


function of centre electrode potential (Va—Ve)?/Va", as measured for 
lenses Nos. 13 and 14, and as calculated by Ramberg for lens No. 14. 


Figure 12 also gives the measurements on lens No. 13 plotted in the same way. 
At low refracting powers this lens is appreciably worse than Ramberg’s lens, but 
approaches this closely at high refracting powers. 
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The influence of the 7 term would be most noticeable for the lower (S— T) 
values ; as the length of a lens does not change with a change of electrode potentials, 
at is reasonable to assume that the values of 7 will not show very great changes, and 
that approximate curves for SD?/f; would be obtained by displacing the 
(S—T)D?/f? curves upwards. 


§7, CONCLUSION 


We may note as a result of this investigation that the main parameter of a 
‘unipotential electron lens is its relative focal length f)/D or refracting power D/fo, 
cand that electrical adjustment of the refractive power has, within certain limits, 
nearly the same effect as an adjustment of the lens dimensions producing the same 
refractive power. Deviations from this rule may be expected for lenses of un- 
usually great Jength, including lenses formed by widely spaced diaphragms, and 
‘lenses which approach the “electron mirror” condition. 

Over the range of lens adjustments which follows the straight-line relationship 
-of Figure 10, the absolute value of the longitudinal spherical aberration Ag* is 
‘proportional to the focal length, and only depends on the value of the constant A, 
and the relative off-axis distance r/D, as (from equation (12)) 


AER, fie ey | eee (14) 


the constant K, being of the order of 8 to 20. Hence the larger the lens diameter, 
‘for a given focal length, the better the lens. 

Electrostatic lenses of the type investigated have spherical aberration 
-constants S which are approximately three to four times greater than those of corres- 
‘ponding magnetic lenses; this appears to be mainly due to the diverging effect of 
parts of the electrostatic lens field, which deflects the rays further from the axis. 
If the diverging effect is eliminated, the electrostatic lens yields aberration 
constants which are very near to those of a magnetic lens of comparable refracting 
power. 

Entrance and exit stops which Jimit the length of the lens field have a consider- ~ 
-able influence on the lens performance, and lead to an increase in spherical 
aberration, approximately by a factor of 2, even when placed in a weak field region. 

The principal surfaces appear to be paraboloids; this fact has to be taken into 
consideration in measurements of the lens constants, or in using simplified 
‘formulae representing the lens aberrations. 

The region of third-order spherical aberration extends in most cases up to 
7/D=0-2 (see Figures 7 and 8), but there are indications in some cases that the 
spherical aberration may decrease again beyond this value of 7/D. An improve- 
ment with wider aperture has been found in the aberration term controlling the 
zonal change of magnification; hence the construction of improved projector 
lenses for electron microscopes appears possible. 
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+ The larger lens diameter explains therefore, amongst other factors, the superiority of the 
smagnetic lens over electrostatic lenses in cathode-ray tube focusing. 
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Age PE IND IX 
I. The “‘ Knife-Edge”’ Test 


From the geometry of Figure 3 we have for the displacement of the shadow 
image of the straight edge 
l—-q, I-q 
Nese pies | atk ae 
foe oe 

with by=q9+ 3d, and b, = Q,Rj ~(q + $d) — (Ag—4Ad) +-1?/2q9. This leads, for 
‘the screen distance ],, to 

__ rylAguld+ $do) — (b— ga) Ad + 472/24) 

* (Go + 34o)(Go + 34 — Ag, + $d, + 72/2qo) 
Further, the magnification M; at the screen in the position J, is 


My=(lj—9)(Qot+ 340). eee (16) 


Hence for 7=1, 2, i.e. two screen positions of which at least one fulfils the 
condition 1;>q), we have four independent equations for the unknowns 
Qo, Aq, $dy, Ad, from which fj, S and T can be determined. 

For /> qo, equation (15) yields Ah, ~7,,I[Aq,, —(4Ad,, +. 72/2q9)]/b2, which gives, 
with Aq, =72S/fy and 4Ad, +72/2q)=72T |fo, for Gyo =hy : 


PR Sr PT Srcas1a) sks bane ate ecTiees a a (17) 


Figure 13. Diagram of distortion figure Figure 14. Optical path through system of two 
obtained in “‘ knife-edge ”’ test. separated lenses. 

From Figure 13 we have, taking the centre of the undistorted image of the 
straight edge as origin A for our coordinates x and y and neglecting the term 
Abie in (Ahh), Ak hxjhy and y?=/=h2-~ But from Figure 3 
h=r(l—q)/fo. Eliminating Ah and h, we obtain the equation of a parabola : 

ten eee 

SS iS SL 

The apex of this parabola is displaced from the origin A by the amount 

Xo =hel (S—T)/(l—qo)?. For y=ho we obtain x=2xy; this is the basis of the rule 

for determining x, and therefore (S— 7) by measuring across the ‘‘ waist”’ of the 
distorted image and along the side of the inscribed square (Figure 5). 

The radius p of the osculating circle at the apex is one half the coefficient of 
in the equation for y*; hence 


(EG Tab > hi 
(S—T)= Bop - (4) Awsewtimvcyest tees (19) 


Il. Addition of Spherical Aberrations of Two Lenses 


Consider two lenses of focal lengths f,; and f, and spherical aberration constants 
SS, and S,, a distance « apart (Figure 14). The focal length of the combination is 
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then given by f=/f,f2/0, where 0=f,+f.—-¢«. From S=Ag*fj/r? =Agf/7? we have: 
r=(Agqf/S)?. From Figure 14 we see that 
ry Aq ) rp (Azle) 1 
tang=>=(s5}) = os f 
eae (Sif f-6 VS: } fa- 

and Age/Agqi =(f2—9)*SolfifoSr. we tees (20) 
The spherical aberration constant S of the combined lens system is then 
S = Agf/r? =AgfS;/Aq fi. But the longitudinal spherical aberration Ag is the 
difference between paraxial focal length f and the focal distance of the ray of 
incidence height 7, ; hence 


Ay pa AeA hae) 


IJ 2 


0 6+ Aq, +Aqe 


For rays not too far from the axis Ag, <0; Ag,<@ hence Ag=(Aq, fo + Agof;)/9- 
Substituting the constants S in place of the Aq’s, and using equation (20), we find 


Sai fs5,--leafaroaee 8) ene (21) 
If one lens is a very weak lens, f, > fo, the expression simplifies to 
S = (Falf) Sate So al Se A eee (22) 


Hence the aberration of the combined system is always greater than that of the 
individual lenses unless the signs of S, and S, are different, a case not normally 
encountered in electron optics. 

Although the constant S, is multiplied by the small factor (/2/f,)’, its influence 
is by no means negligible, since S, increases rapidly with f. If one were to assume 
that the law of Figure 10, S=K,(f/D)?, held over a wide enough range of refractive 
powers to include very weak refracting powers, one would find, by substituting 
this law in equation (22), S=S,[1+(D,/D,)?]. This would indicate that the 
aberrations might be essentially determined by the weak lens as D,>D,. ‘This 
probably overestimates the deteriorating influence of the weak lens, as in the 
limit f, > « the effect of this lens and its aberrations is removed altogether. Never- 
theless, we should still expect an appreciable contribution by the weak lens to the ~ 
spherical aberration of the lens combination even when its contribution to the 
refracting power is quite negligible. 
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ABSTRACT. The change on “aging” in the electrical properties of cold worked 
iron-carbon alloys of low concentration is investigated theoretically. During the process. 
of “‘ aging’, the carbon atoms diffuse into the dislocations, thus forming rows of atoms. 
The electrical resistance due to scattering of electrons by these rows of carbon atoms is. 
expressed in terms of the scattering power of a single.carbon atom dissolved in the iron 
lattice, and hence in terms of the resistance of the same number of randomly distributed 
carbon atoms. It is found that this change depends on the wavelength of the electrons. 
at the Fermi surface. This wavelength is determined by the effective number of free 
electrons, and with a reasonable value of about 0-6 electrons per atom for iron, the agreement 
with the observed diminution in resistance is satisfactory. 

Alternatively, one can use the calculations to determine both the number of free 
electrons in iron and the number of dislocations in the specimen at different degrees of 
hardening, from the observed changes in the electrical resistance and thermoelectric power: 
of the specimen. 


See INDRODUCTION 


T is known that cold worked iron of low carbon content changes its mechanical. 

and physical properties on “aging” at room temperature. For example, 

the critical shear strength for the material increases with time while its. 
electrical resistance decreases. It has been suggested by Cottrell and Bilby 
(1949) that carbon atoms diffuse into regions of the lattice where it is slightly 
expanded by the presence of a dislocation, so as to make the elastic strain energy 
of the dislocation a minimum. ‘They estimate from this the number of carbon 
atoms migrating into a dislocation as roughly one to two per atomic plane. ‘This. 
means (see §2) that the carbon atoms tend to rearrange themselves into rows 
parallel to the lines of dislocations.t It is the purpose of this paper to investigate 
theoretically the change in electrical resistance and thermoelectric power of these: 
alloys, due to such a rearrangement of carbon atoms in them. 


$2, PROBABLE DISTRIBUTION OF CARBON 
ATOMS IN DISLOCATIONS 

The model of a dislocation used by Cottrell and Bilby in estimating the number 
of carbon atoms diffusing into it is shown in Figure 1. In the figure, the «x-axis 
is the direction of glide, and the dislocation is assumed to extend along the z axis 
(perpendicular to the plane of the paper) throughout the entire length of the crystal. 
The slip plane according to the figure is the x—z plane. It should be noticed 
that each row above the x axis contains one more atom than the row below it. 

Using the fact that the strain due to such a dislocation varies inversely with 
distance from the centre, Cottrell and Bilby deduce the extra density of carbon 

* Now at Theoretical Physics Department; University of Liverpool. 


+ The author’s thanks are due to Dr. A. H. Cottrell for supplying the results described in 
§§ 1 and 2 in a private communication. 


PROC. PHYS. SOC. LXII, 4——B 16 


230 A. B. Bhatia 


atoms at any point (7, «) in the lower half of the plane as K sin«/r, where K is 
a constant and « the angle which the radius vector r makes with the slip direction. 
The probability of finding carbon atoms in the upper half of the plane will no 
doubt be extremely small as the lattice in this region is in a state of compression. 
Since the elastic theory breaks down very near the centre of the dislocation, the 
above formula does so also. Nevertheless, it gives some valuable information 
regarding the probable distribution of carbon atoms. Since the probability 
that a carbon atom goes very near the centre of the dislocation in the region A 
(Figure 1) is greatest, it is reasonable to assume that this position will always 
be occupied by a carbon atom in every atomic plane. Now Cottrell and Bilby’s 
calculations show that only about two carbon atoms can go within an annular 
region bounded by circles of radii 2. and 20a. of each atomic plane; and since 
in any atomic plane there are several equivalent positions for the second carbon 


in 


° ° ° ° ° ° ° 


Figure 1. 


atom, one may assume that the positions of the second carbon atom in different 
planes are random. ‘This means that during the aging process, near the edge 
of each dislocation, a single row of carbon atoms parallel to the edge is formed. 
There will naturally be some deviations in the positions of the carbon atoms 
from the straight row. The root mean square of these deviations from the 
straight line, however, will be small in comparison with the wavelength of the 
electrons, whose scattering is responsible for the electrical resistance, and hence 
to a first approximation these deviations may be neglected in calculating the 
probability for scattering of these electrons (see § 3(a)). 

Since a carbon atom occupying an interstitial position in the unstrained 
iron lattice has a tendency to expand the lattice locally, migration of carbon atoms 
into the dilated parts of the lattice due to a dislocation makes the strain and hence 
the dislocation more stable by relieving the stress, but does not appreciably affect 
the magnitude of the strain. Hence any change in the electrical properties of 


the specimen with aging must be due to the rearrangement of carbon atoms in 
the form of rows alone. 
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§3. THE CALCULATION OF ELECTRICAL RESISTANCE 
The calculation of electrical resistance of a metal or an alloy consists of three 
steps: (a) the calculation of the probability that an electron makes a transition 
from a state k to an element dS’ of the surface in k-space having the same energy— 
this is denoted by P(kk’)dS’; (b) the calculation of relaxation time; (c) the 
calculation of. electrical conductivity. 


(a) Calculation of the Probability P(kk’)dS’ 
This probability per unit time per unit volume is given by 


ere 1 dS’ : 
P(kk \dS = 4h . dE|dk, | SeEU Ly dr|?, Sir ete eke (1) 
(see Mott and Jones 1936, p. 251, formulae (14) and (15)), in which %, and 
#,, are the wave functions of the incident and the scattered electrons respectively 
and U is the deviation of the potential from the periodic potential of the ideal 
lattice. ‘The integration is over the whole volume of the material. 

When a metal contains only a few foreign atoms—as in our problem—the 
contribution to U comes only from these foreign atoms; this may be taken to be 
finite only in the sphere whose volume is of the order of the atomic volume. 
Hence denoting by U,, the potential round the mth foreign atom, the perturbing 
potential U can be written as 

Ch ZAG set Lilet Uther 8 rage WTF: eae (2) 


where the summation is over all the foreign atoms. When the foreign atoms are 
distributed at random in the lattice, the value of the integral in (1) is usually 
evaluated at each foreign atom and the total probability of scattering due to all 
the foreign atoms obtained by multiplying by the number of foreign atoms. 
This actually amounts to adding up the intensities of the electron waves scattered 
from different atoms. However, in our problem the carbon atoms in the dislo- 
cations are not distributed at random in the iron lattice and hence we have to 
add up the amplitudes of the electron waves scattered from these carbon atoms. 
For this purpose the integral in (1) may be written as 


{EULY dr = Spi Crib aire nl Td tts oe (3) 


vk.r ik’ 


Now by atheorem of Bloch, },, and 4, are of the forme™"1%,(x, y, 2) ande™ "u(x, y, 2) 
where u(x, y, z) and u,.(x, y, 2) are periodic with the period of the lattice and may 
be taken to be spherically symmetrical for a cubic metal.* We can, therefore, 
replace the integral over the whole volume by a summation over all the carbon 
atoms and integration over a:single atomic sphere. Doing so, we have 
f COP (x U,,)u;,dr = ye *) nl aX, Die dr’ 
nv n 

where dz’ refers to integration over the atomic sphere round the mth carbon atom. 
The integral obviously is independent of the particular carbon atom chosen and 
may be taken before the summation sign. Hence (1) may be written as 


P(kk’)dS’ =fZdS'|Bel™—" me, (4) 


* The magnitude of the wave-vector k in iron may be determined by applying the free ele>tron 
formula to the electrons in the 4s band. The holes in the 3d band are assumed not to contribute 
appreciably to the conductivity since their effective mass is very large in comparison with those 
of 4s electrons. 
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where f2 may be termed the “scattering power”’ of a single carbon atom dissolved 
in the iron lattice as a function of 9, and the summation in (4) is over all carbon atoms. 


f¢ is defined by 
ol adie 
Io | 42?h \dk,, 


We shall assume in what follows that f? does not change appreciably due to strains 
in the iron lattice. 

We now apply formula (4) to calculate P(kk’) for rows of carbon atoms. 
Let the distance between two neighbouring carbon atoms in a row along the 2 axis 
be a. Then if there are m atoms in the row, P(kk’) due to such a row of atoms 
will be given by 


2 
| ut Upuye®—* 8 *n) dr. 


n—1 2 
Py(kk’) = f2|D expfi(k —k’).na} 
0 


=f? sin? {4na(k —k’),}/sin? {Za(k—k’),}, «ss (5) 


where (k —k’), is the zs component of the vector k—k’. Ifthere are N dislocations 
in allin the alloy, with their edges parallel to the z axis and their centres distributed 
randomly in the x—y plane, then P,(kk’) due to N rows of carbon atoms will be 
given by 

P,(kk’) = Nf? sin? {4na(k —k’),}/sin? {ga(k—k’),}. ...... (6) 


In addition the electrons will be scattered by thermal motions, internal strains 
and by any remaining carbon atoms which are distributed randomly throughout 
the crystal. If the number of such carbon atoms be 7, then the probability of 
scattering due to them will be given by 


ft] Tara ie ee eee I oh (7) 


Denoting the probability of scattering due to thermal motions and internal 
strains by P,(kk’), the total probability P(kk’) is given by 


P(kk’) = P,(kkk’) + P,(kk’) + P,(kk’). peered? 


(6) Calculation of Relaxation Time 


The relaxation time for the electrons is given by the well known integral 
equation 


—1=[P(kk (keke dS’, anes (9) 


where the integration over dS’ refers to integration over all portions of the Fermi 
surface, t, and 7, are the relaxation times in directions k and k’, and k, and 
k,, denote the components of k’ and k respectively along the direction of the applied 
electric field F. Assuming the Fermi surface in k-space to be spherical in the 
absence of the electric field, we may to a first approximation write in (9) dS’ = k2dw’, 
where dw’ is an element of solid angle and & is the radius of the Fermi sphere.* 

Let F make an angle « with the z axis and let the plane containing the z axis 
and f be the z-x plane. Remembering that both P,(kk’) and P,(kk’) are 
functions of the angle of scattering alone, it is obvious from the form of P,(kk’) 
that the assumption in the previous sentence does not involve any loss of generality. 

* In addition to the s-s transitions considered abpve, there will be scattering of electrons due 
to s-d transitions. Since neither the mechanism of the s—-d transitions nor their contribution 


to electrical resistance is clear (Seitz 1940), we shall neglect the effect of these transitions on scattering 
of electrons by carbon atoms. ; is 
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Further let k make an angle 8 with the axis and let the plane containing k and 
the z axis make an angle js with the plane containing the zx and x axes and let « 
and ¢ denote the corresponding angles for k’, Then the integral equation (9) 
can be written as 


O18! PRR ge | as Sin p sin a cosd+cosp cos «% i 
nude] ( ue bal ae eae (ee ee at 


BA Re (10) 
In terms of these angles 
dw’ =sinada dd; cos#@=sin« sinB cos(A—y)+cos«cosB ......(11) 
and 
P,(kk’) = Nf; sin? {4nka(cos B—cos «)}/sin? {4ka(cos B—cosa)}. ...... (12) 


In order to solve (10) we must know the variations of f? and P,(kk’) with the 
angle of scattering 0. Since this is difficult to obtain we shall assume that both 
are isotropic, i.e. independent of the directions of k and k’ separately and of the 
angle @ between them. As the potential round any carbon atom may be regarded 
as spherical, and as the contribution to P,(kk’) due to strains is small, the first - 
assumption is justified. Further, since the number of electrons in the 4s band 
of iron is small (~0-6 per atom, see §4) the wavelength of the Fermi electrons 
is sufficiently large to justify the second assumption to a first approximation. 
"Then if we write 


| piinye le pte Pe (13) 
and PPLE ds Hye TS (14) 
we have yok? | Pi(kk’)do!=mq nee (15) 


J 


and 


x1 =F*| Py (ek) deo’ = 4» | 


7 sin? {inka(cosB—cos«)} . 
oO Sg ee 16 

» sin? {tka(cos B — cos «)} ies 2) 

The integrand in (16) has a very strong maximum at «= +f and the limits of 


integration are such that the second maximum at 

TROCCOS B= COS c) = my ieee Se Tt Gs Oe eye (17) 
is not reached since the wavelength is larger than twice the interatomic distance. 
Hence (16) reduces to 


Y=R | Pidsk’ ) deo! =Nan(n/hka)=nyy(n/ka), ee . (18) 


where n, denotes the total number of carbon atoms in all the dislocations. It 
is seen that none of these integrated values depend on the incident direction k. 
An occasional gap in a row of carbon atoms will not affect the results, since so 
long as the length of the rows is large compared with the wavelength, formula 
(18) holds and x, depends only on the total number of carbon atoms that have 
gone into the rows. 

The integral equation (10) can now be written 


1 AR? j : ry fsinp sina cos d + Cos pu Cos « ; 
Bears a x | {B(ie!) + Pa(leke’) + Pa(R’) yr ‘an usin B cos + cos pu 4 


where x= x1 + X2+X3- 
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This is the Fredholms integral equation of the second kind and its solution 
can be obtained by the Liouville-Neumann method in the form of an infinite 
series (Margenau and Murphy 1943). The evaluation of different terms in the 
series becomes particularly simple in view of the assumed isotropy of P,(kk’) 
and P,(kk’). ‘The infinite series obtained as a solution of (19) is convergent 
and can be summed giving 


poe a 
n= = [14% sin p sin 8 cos s+ cos p cos B Litt x e x aaa 
=<[1+ A oa eo | Bas os (20) 


X—X1 sinp sin 8 cosy +cos pu cos B 
since y;<x. ‘The solution may be verified by substitution in equation (19). 


(c) Calculation of Electrical Conductivity 
The specific conductivity is given by (Mott and Jones 1936, p. 262, equation 


(52)). aa oe 
a Jbsaalgel on GAZE 
he an | | a) Tp teadh Ak. yes 8) 

Taking the spherical polar coordinates for k, as before, (21) becomes 


oe (ie Td Ones sors. 
o= — Fe aap | | |r gy (Sin sin B cosy+ cos cos f) an sin B dB dis dk. 


Since (dfo/dk) is finite only in a very small range of k at the surface of the Fermi 
distribution, while integrating with respect to k we keep all other terms as constant 
at k=k,,,, and, with f(df)/dk)dk = —1, we have 


e dE ; 
= Ot : : te 
oe ane (: ah iy {| 7 (sin sin 8 cosy +cosp cos B)* sin BdB dbs 


2 rT 2st 
= uae ee | |  -7,(sin pw sin B cos +cos p cos 8)? sin BdB dy, 
M1 An B=0 0 


y= 


So dBi: (22) 
where N.¢ is the effective number of free electrons per unit volume * and is given 
by 

5 1 
New = {5 (=) aS Me Paes (23) 
(Mott and Jones 1936, Chap. III, §4). 
Now substituting for 7, from equation (20) and integrating, we have 
a Nese ; d E ae cos? | a 8 RS TUE TE (24) 
Git aX Xk 


It is obvious from (24) that when .=0, that is when the direction of applied 
electric field is along the z axis, the conductivity is a maximum and equal to 
Nwe*/m(x2+ x3). This means that the rows of carbon atoms do not offer any 
resistance to an electric current flowing along their direction. This can be seen 
from general considerations also, as follows. It is clear from expression (12) 
for P,(kk’) that, for a given direction k of the incident electrons, P,(kk’) is finite 


* Neff will be equal to the number of electrons in the 4s band of iron. 
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only in the neighbourhood of those directions of k’ for which the z component 
of the momentum remains unaltered. Hence if the electric field is applied 
along the z direction the scattering by the rows of carbon atoms will not obstruct 
the flow of the current in this direction. 

It can also be seen from (24) that the conductivity is least for ~=47, when 
o=N.,e?/m x or the specific resistance 


R=(m/N pre”) x = (m/Nesre)(X1 + X2 + Xs) 
is a maximum. If all the carbon atoms present in iron go into rows then for 


b= dT, 


R= (m|Nee2)(x1 + x3); 


as compared with R=(m/N.e?)(x2+ x3) when the carbon atoms are distributed 
at random. Now for the same number of carbon atoms 


(vas) ae Ces (at eee aa ck (25) 
and since N.g~0:5 per atom (see §4) and k can be obtained from the formula 
Sie = (2a) New ae ree (26) 


we havec~1-3. This would correspond to an increase in the resistance for current 
in this direction. 

On the other hand if one is experimenting with a polycrystal or a very much 
hardened crystal, so that the rows of carbon atoms can have all possible orientations 
with the direction of the applied electric field, then the observed electrical resist- 
ance will be the average given by 


(fe | we x(x=%) ee] sin du de 


R= _m_Jolo Lemma conte | ne de Xa 4X OS" ee {(2 — 1) tan ( — ie é 
av nS N. ee 
ef ee if in sin p du de eff X1 X1 
0/0 


§4. COMPARISON WITH EXPERIMENTAL -DATA 
In a cold worked material, there are about 10" or 10!” dislocation lines per 
unit area (Seitz 1943). Since, as we have seen, one can assume that on aging 
only one row of carbon atoms is formed in each dislocation, the number of carbon 
atoms taking part in the formation of rows will be given by the following table: 


No. of dislocation lines per unit area nie Oz Ose 
No. of carbon atoms per unit volume i B43). <0" Sh4 WOES 
No. of carbon atoms in atomic per cent .. 0:46 x 10 4-6 x10? 
No. of carbon atoms in weight per cent .. 0-001 ~0:01 


The number of atoms per cm? in iron are 8-55 x 10°. 

Now 0-01 weight percent carbon increases the electrical resistance of iron at 
room temperature by 2:5°% (Dijkstra 1947). Hence, using the fact that 
X1/(x — x1) <1, we have from equations (24) and (27), 


R=(m/Newe)(Xs+X2+x1 SiN?) wee ee (28) 
: ne ioc GIN ROAVaNtIE, FP Esc Se ame (29) 
and Reet oot ewe OT eee (30) 
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where R, is the specific resistance due to thermal vibrations, internal strains and 
due to any carbon atoms distributed at random, and R, sin? is the resistance 
due to rows of carbon atoms. If R, is the specific resistance of the specimen 
before aging, i.e. when all the carbon atoms in the specimen are distributed at 
random, we have from (25) 

Ro=Ra+ Role eon ae Sige a4 eas (34) 


Tn Figure 2, Ry and R are plotted in arbitrary units against ~. Curves (a) and (6) 
represent the values of R when N.g.~1 per atom (c~1-1) and when Ng ~ 0:6 per 
atom (c~1-3). The origin refers to the value R,. 

For the case in which the directional distribution of the rows of carbon atoms 
is random, we can express the diminution in resistance in the form 


R, =AR(1 —2c/3) PR Ey) 


where c is given by (25) and AR is the resistance due to the same number of carbon 
atoms distributed at random as go into the rows on aging. 


R 


Figure 2. 


In order to be able to compare the value of R, as obtained from (32) with 
experimental values, we have to know c. It is by no means easy to decide a 
value of Nw for iron. Manning (1943) has made some calculations for the 
number of electrons in the 4s band of iron, or Ny, as 0-26 per atom. ‘This value, 
however, seems to be too low, since the binding energy and electrical conductivity 
of iron are of the same order of magnitude as those of nickel and cobalt, for which 
the number of free electrons is ~0-6 per atom (Mott 1935). From the data on 
saturation magnetic moments of alloys of iron with Mn, Cr, Co and Ni, Slater 
(1936) also was led to assume the same value of 0-6 per atom for N. for all the 
three elements Ni, Co and Fe. We shall, therefore, assume that Nog for iron 
is 0-6 electrons per atom. Using (26) we have c~1-3 and R,=0-4AR/3 as 
compared with c~w1-1 and Ry=0-7AR/3 if Nog~1 per atom. 

Assuming that there are 101” dislocation lines per cm?, then since AR =2:5°%, 
of the resistance of pure iron, Ry =0-3% of the resistance of pure iron, as compared 
withthe observed value of about 0-2°%. Thisshould be regarded as a fair agreement 
with theory, in view of the various simplifying assumptions involved in the theory 
and an arbitrary assumption about the number of dislocations. 

It may be pointed out that if one uses a single crystal which is worked in ee 
a way that all the dislocation lines are parallel to each other, one can determine 
both c and the number of carbon atoms gone into the dislocations to form rows, 
- and hence the number N of dislocations themselves, from (29) and the observed 
variation of the specific resistance with the directions of the applied field. It is. 
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hoped that experiments of this type on single crystals when carried out will give 
-some valuable information regarding the number of dislocations at various stages 
-of hardening. 


$5. THERMOELECTRIC POWER 


It is clear from expressions (15) and (18) that the electrical resistance of carbon 

-atoms depends on the wave number k (or the energy £) of the electrons at the Fermi 
surface in a different manner, according as the atoms are arranged in rows or 
are distributed atrandom. Hence, since the thermoelectric power is proportional 
to d(log R)/dE, it is obvious that on aging the specimen will show a change in the 
thermoelectric power. This change naturally will vary with the direction of 
the applied electric field when all the dislocation lines are parallel in the crystal, 
and hence can be utilized to determine c and N, provided one knows the dependence 
-on R (or E) of the resistance of pure iron and that of the carbon atoms distributed 
at random in it. The required dependence can be easily obtained from the 
measurements of Dupuy and Portevin (1912) on the thermoelectric power of 
iron of low carbon content. 

For a polycrystal also c and N can be determined by comparing the observed 

changes in both the specific resistance and the thermoelectric power with suitable 
equations from the theory. 
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ABSTRACT. A calculation is made of the change in the energy of a monovalent metal: 
when one of its atoms is replaced by an atom of a different monovalent metal, not necessarily” 
of the same atomic volume. From a knowledge of this energy the heat of solution can be 
determined and thereby, through the usual thermodynamic relations, also the solubility 
limits. Numerical calculations have been carried out for the copper-silver system and the 
theory shows in this case the existence of narrow solubility limits, in agreement with the 
known phase diagram. The theoretical values of the heats of solution are compared with 
those required to give the observed solubility limits. 


yf 


Si INTRODUCTION 


N this paper, an attempt is made to calculate the heat of solution when one- 
monovalent metal is dissolved in another and thus to find, through the usual 
thermodynamic relations, the position of the solubility limits where such limits. 
exist. ea 
In the case when the constituent metals have different valencies it has been 
shown that the Fermi energy, determined by the Brillouin zone structure, is mainly 
responsible for the solubility limits and this led to a qualitative explanation of the: 
Hume-Rothery rules (Jones 1934). However, for the case in which the solute and 
solvent metals have the same valency, very little work on the calculation of the 
energy of an alloy as a function of its composition has hitherto been carried out. _ 
The copper-silver system was chosen for detailed calculations because, in the. 
first place, solubility limits are known to exist, and secondly because Hartree self-- 
consistent fields of the copper and silver ions, which are required in the calculations, 
are available. 

The following treatment is based on the general ideas of the Wigner and Seitz. 
theory of the constitution of pure metals. ‘The main difference arises from the fact. 
that the potential in the alloy is no longer periodic and this requires a somewhat. 
different method of calculating the Fermi energy. 

In §2 the heat of solution is defined and a relation between the temperature: 
and the atomic concentration is obtained. In §3 a method for obtaining one- 
electron wave functions in an alloy is described. The change in the Fermi energy 
of a metal when a foreign atom replaces one of its atoms is dealt with in §4 and §5. 
The change in the total energy is considered in § 6. 


§2, THE HERAT OF SOLUTION 


The energies of all the metals considered are defined with reference to the same- 
zero level of energy, which is taken to be the energy of all the ions and electrons at 
infinite separation. 

Consider an alloy of two metals A and B. Let c, and c, be the atomic concen- 
tration of B in two phases of the alloy, and f, and f, be the free energies per atom of 
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_ these phases respectively. The conditions for the equilibrium of these phases, at a 
. common temperature T, are 0f,/dc, =0f,/0c,=(f, —f2)/(C, Ce). If c, denotes the 
_ heat capacity per atom at constant volume, and E(c) the energy per atom at the 
_ absolute zero, the free energy per atom of a phase at concentration c and temper- 
< _ ature T may be assumed to be given by the following eee 3 


fle, T)=hT{elne+(1—e)ln(1—0)} + (0) + | c,dT — T| aes aT/T. 


| Applying the previous equilibrium conditions and assuming thas Hie change inc, 
is negligible when c varies, it can be deduced for dilute solutions that the concen- 
tration c of metal B is given approximately by 


kT Inc=E(1)—E(O)—[@E(c)/Ac]..g eevee (1) 


D where E(1) and E(0) are the energies per atom of the pure metals B and A respec- 
tively. 

If Ng denotes the number of atoms of type B and N the total number in a given 
lattice, then dE dc=[6(NE) ONy]xy and, therefore, [GE/dc], ,, denotes the change 
in the energy of the lattice, of a dilute solid solution of N atoms, where one atom of 
_ A is replaced by one atom of B. Hence W=E(1)—E(0)—[0E/dc],_, is the 
__ energy liberated when an atom of metal B is taken from the pure metal B and put 
| into the dilute solid solution of Bin A. The quantity so defined is called the heat 
of solution. Phase boundaries exist when W is negative, that is to say when work 
must be performed to add a solute atom to the solid solution. The quantity 
E(1)— E(0) is known, since for either metal E= —(S+J) where S is the heat of 
sublimation (Grimm and Wolff 1934) and J is the ionization potential of the atom; 
hence we can determine the solubility limits if we can calculate [0E/dc],. 
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$3. THE VALENCE ELECTRON WAVE FUNCTIONS 
Consider a lattice of copper atoms in which a single copper atom is replaced by a 
silver atom. Both copper and silver have a face-centred cubic structure and it will 
_ be assumed that the interatomic distance d of the twelve nearest neighbours of an 
atom is not changed by the presence of the foreign atom. Ifata point ata distance 
* from the silver atom, along each line of the nearest neighbours, a plane is drawn 
perpendicular to that line, then these twelve planes will enclose a dodecahedron 
around the silver atom. Let the position vector of the silver atom be ry and that 
' of a copper atom be r,. Inside this dodecahedron the potential energy of an 
electron is taken to be V(r — 14) + C(r) and outside it, in the neighbourhood of the 
copper ion at r,,, which is nearest to r, the potential energy is taken to be V,.,(r—r,). 
V,, and V,,, are the potential energies of an electron in the fields of a single silver 
ion and a single copper ion respectively, defined so that they vanish as 1/r at large 
distances, and C(r) depends on the distribution of the electronic charge round the 
silverion. If C(r) is known, then for a given choice of the distance x the potential 
energy of an electron anywhere in the lattice is determined, and may be denoted by 
V(r). 

If the one-electron wave functions could be found in the field of this potential, 
then the total wave function could be written as a product of two determinants, and 
the expectation value of the exact Hamiltonian operator calculated. The 
Hamiltonian operator H is ae by : 


H= z Hilt) BE z ee hone (2) 


545 °F a>bl ab 
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where H,(r;) = —V?7+XV,(r,), x; is the position vector of the 7th electron, 7;; the 


a 
distance between the 7th and jth electrons, 7,,, the distance between the ions a and 
b and N is the total number of electrons in the lattice. In this expression the 
unit of energy is taken to be the Rydberg unit, viz. 13-53 ev., and the unit of length 
the radius of the first Bohr orbit. 

The best values to take for x and C(r) would be, according to the variation 
principle. those which minimize the total energy calculated in this way. As a 
zero approximation C(r) has been assumed to be equal to zero, and the distance x 
to be fixed by the condition 


Vi (SY Pat ce ek ee eee (3) 


It will be shown that there is either an accumulation or a diminution of charge 
round a foreign atom and thus a value of C(r) which is different from zero would 
lower the total energy relative to the zero approximation. If ¢g, denotes the excess 
charge at the foreign atom, then C will be of the order of magnitude q,/7) where 7p is 
the radius of the atomic sphere. In some cases, this will be far from negligible, 
and it would then be necessary to proceed to higher approximations. 

The one-electron wave functions % with energy Z, have been obtained in the 
following way. Write Y(r)=y(r). u(r) and E=e+e 9, then the wave equation, 
expressed in Rydberg units of energy and lengths in units of the first Bohr orbit, 
can be written 


ulV2x + ex] + xf V2u 4 fey — Vir) a] = —2V x. Vu. a. (4) 


For a pure metal V(r) is periodic, and an approximate solution can be obtained, 
first by neglecting the term 2Vy. Vu, and secondly by equating both brackets on the 
left hand side to zero. ‘This gives a solution of the form %=u(r).exp(k.r), 
where k? =< is the energy of an electron relative to ¢) which is the lowest eigenvalue 
corresponding to a solution in which u is periodic with the period of the lattice. 
With this approximation the total energy is given as that of the lowest state ey plus 
a Fermi energy given by the free electron approximation. 
An attempt is made now to determine the wave function of an electron in the 
copper lattice in which a silver atom replaces one copper atom. Defining a 
function u(r) such that v(r) =¢ inside the polyhedron round the silver ion, where 
‘p is a constant, and v(r) =0 outside this polyhedron, equation (4) can be written 


u[V2x + {ce —v(r)}v] + x[V2u + {eg + u(r) — V(r)}u] = —2Vy.Vu. ...... (5) 


‘The Wigner and Seitz treatment of pure metals suggests that, as a first approxi- 
mation, the term 2Vx. Vu should be neglected and solutions of the following two 
‘equations sought : 


Vea feg PUT) — V(t) a ly Pee (6) 
and Veveies Of) | peed o's matomba acess (7) 


‘Taking for the domain of the wave function a large sphere of radius Z in which 
there is a single foreign atom at the centre, the boundary conditions are : (i) 
a vanishes over the sphere of radius L, (ii) # and #’ are continuous everywhere 
inside this sphere. 

To get the solutions of (6) outside the central polyhedron, that is where v(r) =0, 
each polyhedron round a copper ion is replaced by a sphere of the same volume and 
the boundary condition u’(r))=0 is applied over this sphere of radius ry. This 
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approximation of replacing the polyhedron by asphere leads to asmall discontinuity 
in the derivative of u(r) at the surfaces of the spheres, but Wigner and Seitz have 
shown that the error involved is negligible. 

The solution of (6) inside the central polyhedron is obtained by replacing this 
polyhedron by a sphere of radius a which has the same volume, and using the 
conditions of continuity of u and w’ at the surface of this sphere. It is clear from 


similarity that 
a|x=19)$d = (2r,—a)|(d 2), 


and since the potential of the ion at distances of order a or 79 is practically the 
Coulomb potential it follows that 


V s.@)/V cal279 — 4) = (219 — @)/a =(d — x) /* = Vig g(x)/Vey(d — 2). 
Thus condition (3) can be replaced by 
Va hee Gall or Grr setts uy a Tee Oe (8) 


‘Taking the values of these potentials to be those given by the Hartree field (Hartree 
1933, Black 1934-1935), (8) determines the radius of the central atomic sphere a. 
Within this sphere V(r) = V,,(r) and u(r) =¢ which is determined so as to assure 


5 -0-01 
At ey nl 7; Nl aS SS ae ae ees) 
24 28 32 36 2 35 31 27 23 19 
T measured from solvent atom r measured from solute atom OU 
Figure 1. Full line represents the radial wave Figure 2. Illustrating the method of 
function in the alloy of Cuin Ag. Dotted applying the boundary conditions 
line represents the radial wave function in for the determination of the wave 
the Ag pure metal. function of the lowest state. Curve 


(a) 6=—0:75. Curves (6), (c), (d) 
¢=0-05, 0:063, 0:08, respectively. 


the continuity of the wave function and its first derivative at the surface. Solving 
(6) for different values of ¢ a family of curves of [w’(r)/u(r)],, can be plotted 
and the solution which satisfies the condition 


u'(a) u' (279 — a) 
ae oe io sees (9) 
u (a) Ag Uu ( Yo a) Cu 
will be continuous to the same kind of approximation as the wave function 
obtained by Wigner and Seitz method for pure metals. Figure 1 shows the form 
of the wave functions obtained and Figure 2 shows how ¢ is obtained from the 
condition (9). It is found, in the case of a silver atom in a copper lattice, that 


1) = 3°3, a=3-34 Bohr units and d =0-063 Rydberg units. In the case of a copper 
atom in a silver lattice 7) =3-75, a=3-72 Bohr units and ¢ = — 0-046 Rydberg units. 
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Equation (7) is the Schrédinger equation appropriate to the problem of a 
particle in a sphere of radius L in which the potential is zero everywhere except 
over a sphere of radius a, where it is constant and equal to ¢. The solution can be 
written y(7)=r ?R(r)Y_,,,(8,¢), where 

R(r) = AJ s(7V (e— 4)), (0<r<a), 
and R(r) = BI (1V/e) + CI (rv); (a<rcL). 
R(r) in the range (0 <r <a) is denoted by R,(r) and in the range (a<r<L) by R,(7). 
The conditions R,(a)=R,(a), R{(a)=Ri(a) and R,(L)=0 determine the eigen- 
value e. Let N,(e)de be the number of eigenvalues in the range from « to « + de for 
given J, then the total energy of the electrons in the lattice relative to that of an 
electron in the lowest state will be 


Ey =d(21+1) |” N((e)ede, 


where e, is the maximum energy of an electron in the Fermi distribution relative 
to that of the lowest state. £), is referred to as the Fermi energy. 


S45 TEE GLA NGE SING DAES bERM lb NERG YS Di@E eos EL, 
PRESENCE OF A FOREIGN ATOM 


If < is the eigenvalue in the solution of (7) for the pure metal, in which 
v(r)=0, and E{ the corresponding Fermi energy, then the change in the Fermi 
energy when a single foreign atom replaces an atom of the pure metal is given by : 


By ~ Bi) =¥2(21+1) ie N(e)[e —< de. 


‘This assumes that the density of states is not affected by the potential ¢ in the 
central sphere of radius a. For /=0 the eigenvalues are given by the equation, 


_ ff tanavV/e— vV/{e/(e—¢)} tan av/(e—¢) 
tan. Lv/e as Fepacesncrevaeseel Oe” Vel ‘ater eters (10) 
When ¢ is positive and larger than e, tan[ay/(e —4)]/+/(e —$)i is replaced in this 


equation by tanh ae ei mee) When ¢=0, i.e. in the case of a pure 
metal, the eigenvalues «© are the solutions of the equation tan (L/e) =0, which are 


EO) = 22] LE? BES eae” FL eee (11) 
Since L is very large the left hand side of (10) fluctuates rapidly as « increases 
whilst the right hand side is a slowly varying function of «. Hence, without 
involving any essential approximation, « can be replaced by < in the right side of 
(10). If the right hand side of (10) is denoted by tan {(<), then 
Lyfe = bo(e) + n7, i §t lected 
and using (11) the change in the eigenvalue of (7) may be written : 
€— €O = 27-3 4/ML (€). meote eldacd)) 


Equations similar to (10) are also obtained for values of / different from zero; for 
J=1 the equation is 


— glacotay’« + («—~o)vVe tev (e—$) cot aye cotay/(e~$) 


anelia/e= 
tan a/ p/a+(e—¢)v/e cot av/e—€/(e— 9) Cot av/(e— 4) 


ee oe el 


P, ea oir Lacie 
A A tl ye 


ee 
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In each case let tan Ze) denote the right hand side of these equations. Then 
tan Ly/e=tan(e) and Ly/e=G(e)+nz7, and the total change in the Fermi 
energy can be written 


Ey—E® =32(21+1) i BN OiRigfeleyier | ha (13) 
1 0 

In this summation the contribution to E,—£\° of the terms corresponding to 

small values of / may be written 


¥2(21+1) ie Sh 0 ee eee 
1 0 7 


since for small values of / there is no appreciable difference between N,(e) and N,(e) 
which is equal to L/27./e. By calculating numerically the sum (14), it is found 
that the terms corresponding to /=0, 1 and 2 are the most important, the others 
being negligible. Figure 3 shows the curves of the functions ¢,(e) for these values 


: € 
0:0435 0-0870 0:1305 0:1740 0:2175 0261 
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Figure 3. Displacement of the energy levels Figure 4. Correction to the eigenvalues given 
of the conduction electrons as given by by (19). Curve (a) refers to the silver pure 
equation (11a) due to a copper atom in metal and (4) to the solution of copper in 
solution in a silver lattice. silver. 


of 7. For one silver atom dissolved in a copper lattice it is found that the change 
E,—E is 0:0588 Rydberg units, and when a copper atom replaces one of the 
atoms of a silver lattice the change is —0-0484 Rydberg units. 


§5. CORRECTIONS ARISING FROM THE TERM 2VU.Vx 


In the previous section the eigenvalues of equation (4) have been obtained by 
neglecting the term —2Vy.Vu. The effect which this term has on the eigen- 
values, and therefore on the total energy of the electrons, may be estimated by an 
application of the perturbation theory. If Ae denotes the increase in the eigen- 
value due to the term —2Vy . Vu then, to the first order, 


Ae = = | xu. 2x. Vide | | ude, seen (15) 


This quantity will be shown to be of order 1/Z and therefore makes no significant 
change to the energy of a pure metal. As Wigner and Seitz have shown, for the 
pure metal, it is necessary to consider the second order correction to the eigen- 
values. Here, however, it is the change in the eigenvalues, brought about by the 
replacement of a solvent by a solute atom, which is of interest. ‘This change is 
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itself of order 1/Z. Thus, in this case, even the first order term (15) will lead to a. 
significant correction to the difference of the eigenvalues just mentioned. 
In order to reduce the integrals of (15) first let 


1 
u*(r)='—— || wrydSy enews 16 
= pal woas, (16) 
where the integration is over the surface of asphere of radiusr. This function of r 
must tend rapidly to a constant value, c?, asyincreases. In fact it may be assumed 
that wv? will only differ appreciably from c? within the central atomic sphere. The 
denominator of (15) may now be written 


‘ ih SA cr 

| eu dr =4r | Pyrdr sana (17) 
0 

It has been shown that the main contribution to the required energy change comes. 

largely from those states for which /=0. The corrections required by the term 

2Vx.Vu will, therefore, be calculated for these states only. For these states 

x =(1/r) sinr+/e, and (17) becomes 


-L <= 
i ud =4n | (sin? r/c 2(r) dr, 
0 


and since it is only in the first atomic sphere that uv? is at all different from c?, this 
integral is equal to 2mc?L without any significant approximation. 

The integrand of the numerator of (15) may be written (yVx)Vu? and applying 
Green’s theorem (15) becomes 


1 
he= = | VGN) ne (18) 


This integral is taken over the sphere of radius L, and using (16), (18) reduces to: 


new 2 |’ PO) 4 (yf 4%) 
dx 


L 
Since J, 5 (x Nee o) dr=0, (18) may be written : 


hers ane -1}2 (7 x2) dite Valea (19) 


The integrand of (19) becomes vanishingly small when 7 is appreciably greater than 
the radius of the central atomic sphere. It is assumed, therefore, that Ae may be: 
estimated by replacing the upper limit L by ain the integral of (19). When wand 
x refer to the wave functions for the pure metal the energy correction is written 
Ae. When a solute atom occupies the central atomic sphere the appropriate 
functions u and y, as calculated in § 3, are used in (19) to give Ae, and the correction 
to the energy difference is thus Ae — Ae. Figure 4 shows Ae and Ae calculated 
in this way for pure silver and for a silver lattice containing a copper atom respec- 
tively. The difference Ac —Ac® will be written 2(¢)/L. In this calculation 
no(e) has been determined only for the states /=0. For any other / the corre- 
sponding difference may be denoted by 7,(e). ‘The corrected change in the Fermi 
energy will be (E,— EF) +AE,, where 


AB y =S2( 21+1) ["* =1) Eee = 
é 
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| For a silver atom in a copper lattice the value of AE, is 0-0063 Ryd., and for a 
| copper atom in a silver lattice AE, is —0-0133 Ryd. 


$6. THE CHANGE IN THE -LOTAL ENERGY 


As a zero approximation, it has been assumed in §2, that the potential C(r), 
due to the charge distribution round the foreign atom is zero. The calculated 
wave functions, based on this assumption, lead to a considerable change in the 
charge inside the central atomic sphere if the enclosed atom is replaced by a 
foreign one. Ina silver lattice with one copper atom, the wave function outside 
the centra] atomic sphere is of the form 


o(t) = Ar sin {(r — L) v/e} ua,(r); 


and inside this sphere 
-0(7) = Br* sin {7-/(e — $)} Uou(7). 


From the continuity conditions it can be proved that 


A®/B? =[e +4 cost av/(e—$)] |e, 


and the value of the normalization factor A can be taken, as stated before, to be: 
given by 1/A?=2nc?L. Taking the density of state N(c) to be L/27/e, the 
change in the central charge may be written as : 


ued &p 4/e[1—cos 2ra/(e— p) |de a se g 
later |. er ccoete opi ~ =, fVer—sin @rvex) 27] 


(20) gives an accumulation of charge of 32°% round a copper ion in a silver lattice, 
and a diminution of 29°% in the charge rounda silver ion inacopperlattice. This. 
change in the charge distribution leads to a change in the potential energy of the 
lattice. If in a higher approximation C(r) is taken to correspond to the altered 
charge given by (20), a different change in the charge distribution will be obtained. 
However, this second approximation would only slightly affect the change in the 
total energy, since the change in the Fermi energy would be affected in the 
opposite sense. 

To get an expression for the total energy of the lattice, it is sufficient to take the 
total wave function of the electrons as a product of N one-electron wave functions 3. 
since the change in the correlation between electrons is being neglected. The 
total wave function is taken, therefore, as 


V=,,(t1).--- Pc y(Ey)> 
where y, denotes the one-electron wave function satisfying the equation 
Vitac(ti) + [ex — 2 U a(t) Yx(ti) =9, 


and is normalized so that {| #,(r) ?}dr=1, over the sphere of radius L.  U,(r;) is 
defined so that U,(r;)=V,(r;), inside the polyhedron round ion a and is zero 
outside it. Using the Hamiltonian operator given in equation (2), we have 


| we | PDH (H)Y dry. dry = Dey, + EE | PE (EV a(t,)— U alt) (Fier. 
a a € ad 
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The expression =| GPRD mc U,}t;,47; gives the potential energy of the charge 


distribution in the whole lattice except inside the a cell, due to the field of the ion a. 
Outside the a cell V(r) — U,(r) = —2/r,, where 1, is the distance to the ion a, and if 


the lattice is imagined to be divided into cells of equal volume {2,, around each ion a, 


then 


ba 


sa: ide: 
E | de(Va-Une,dri=— % | =ple)ar, 
a4 Q,7 a 


where p(r) is the charge density at the point r. The part of the potential energy 
of the lattice due to the electrostatic interaction between the electrons is given, by 


bole 


2 2 
BS | = plry)olms)drydry+ BE | | ples) (es) drs dre 
a 49,19, 712 a#b/ 2,1 2,712 


The first term is the self energy of the charge distribution, and the second term is 
the mutual interaction energy of the set of polyhedra. Thus the expectation value 
of the total energy in the lattice can be written ~ 


; 2 
<H) = Ey +¥t0,—ZE = elt) dr +B :| {= = elts)p(te) dry dra + EE, 
Q, Q, 


aHxb QT axb a+b Tab 
where the last term arises from the interaction between the ions, Ey = Ze,. is the 


Fermi energy and dw, is the self energy of the charge distribution. Supposing 


a 
that each cell can be replaced by a sphere within which the charge distribution is 
spherically symmetrical and (1 + q,) is the total negative charge inside the a sphere, 
then it follows that 


a = Ey ah LW, as ee = I aw JaIv)|T av: SS 9. sien (21) 


38 a uniform charge distribution is assumed inside each a sphere, w,, will be equal 

to (3/5)2r) 1g where 7, is the radius of the sphere. The change in the total energy, 
when one atom of a pure metal is replaced by a foreign atom, is denoted by 
<H> — <H®>, where <H> corzesponds to a charge distribution with g,=0. 
<H> — <H®> is to be compared with [0E/dc],_) obtained from the experimental 
values of the heat of solution (Ageew, Mansen and Sachs 1930). 

The forms of the wave functions, obtained with the approximations of the 
present calculation, suggest that the excess qg in the charge inside the central atomic 
sphere, is accompanied by a decrease q/N in the charge inside every other sphere. 
As it is difficult to calculate the exact charge distribution, the above distribution will 
be considered as a limiting case. The last term in (21) is then vanishingly small. 
For the solution of silver in copper the calculated value of <H> — <H> is 
0-0982 Rydberg units, while [0Z/dc],_.» =0-0707 Ryd., for the solution of copper in 
silver <H> — <H> = —0-0262 Ryd., while [@E/dc],_, = —0-0324 Ryd. 

Another limiting case, in which the excess in the central charge is compensated 
by a decrease in the charge of the immediate neighbourhood, is suggested by the 
calculation of Huang (1948). ‘The actual distribution of charge must be inter- 
mediate between these two cases. With this second assumption that only the 
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charge round each of the twelve nearest neighbours of the foreign atom is affected, 
it is found by direct calculation of (21) that <H >— <H®> is 0:0817 Ryd., when a 
silver atom is dissolved in a copper metal, and j is —0-0439 Ryd., when a copper 
atom is dissolved in a silver metal. 

The experimental values for the heats-of solution of Ag in Cu and Cu in Ag 
are —7-2and — 6-8 kcal/molerespectively. The corresponding theoretical values, 
according to the above figures, are between —16-4 and —10-4 kcal/mole for the 
first case Ag in Cu and between — 4-8 and —0-3 kcal/mole for the second case Cu in 
Ag. Although there is, therefore, no close agreement between the calculated and 
observed values of the heat of solution W, it should be remembered that W is 
obtained as the difference of two nearly equal quantities with the consequence 
that errors arising from the inevitable assumptions and approximations are greatly 
magnified. ‘The fact that the theory leads to negative values of W and therefore to 
limited solubility in both cases is, in itself, of interest. It would bea searching test 
of the method, to apply it to the copper—gold system where the solid solutions are 
continuous. ‘The reason this calculation has not already been made is that there 
appears to be no sufficiently accurate self-consistent field for gold yet available. 
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Mechanical Breakdown of a Non-Metallic Polycrystalline 
Material Observed by Reflected Light Microscopy * 


By J. P. ROBERTS 


High Temperature Materials Group, Metallurgy Department, 
Royal Aircraft Establishment, Farnborough, Hants. 


Communicated by L. G. Carpenter ; MS. received 20th October 1948 


ABSTRACT. Astudy of fractures in a well-sintered polycrystalline aggregate of corundum 
crystals has shown that the strengths of the crystals and the crystal-boundaries are approxi- 
mately equal. Reflected light microscopy was employed and the usefulness of this method 
of examination has been established. 


Si ENERO De CL VON 
LOOSE-PACKED polycrystalline aggregate of corundum (a-alumina) 
crystals may be formed by slip-casting+t or pressing. The specific 
gravity of such a mass is commonly 1:8—2-5, whilst that of a corundum 
single crystal is 4-0. If the aggregate be maintained for a few hours at 1,650— 
1,800°c. (melting point is 2,050°c.), some crystals grow at the expense of their 
neighbours and the process causes a diminution of the amount of void space and a 
corresponding increase in the specific gravity which may reach finally a value 
within 2°% of that of the single crystal. ‘The process is known as “‘sintering”’. 
The mechanical strength of polycrystalline corundum bodies made in this way 
is a function of the temperature and duration of sintering. Before sintering, the 
aggregate is fragile and possesses negligible strength, but as sintering proceeds the 
strength increases, due to an increase in the strength of the bonding between the 
crystals and a reduction of the amount of void. Reliable information concerning 
the relative strengths of crystals and crystal-boundaries has been lacking. This 
paper describes a microscopical study of fractures in well-sintered corundum 
bodies; knowledge of the relative strengths has been gained. ‘ 
Transmitted light examination was not adopted for observation of fractures,. 
since, firstly, preparation of the necessary very thin sections in a plane perpendicular 
to a fracture is not fully satisfactory, due to the difficulty of avoiding loss of material 
during grinding, and, secondly, crystal-boundaries cannot always be discerned 
easily. Reflected light examination of this class of material has not been practised 
in the past, no doubt partly due to difficulties in polishing; corundum follows 
diamond on Mohs’ Hardness Scale. A technique of polishing and etching 
sintered corundum was worked out for the fracture study. 
Some observations concerning crystal-boundaries and etch-figures are included. 


§2. EXPERIMENTAL PROCEDURE 


Sintered corundum bars with crystal-size of approximately 1004 were 
fractured in bend (stress in outermost fibres at fracture=12-17 ton.in-?, ice. 
1-9-2:6 x 10° dyne.cm). The broken pieces were then fitted together, taking 
care that no parts of the fracture faces were damaged in this operation. After 
embedding in a suitable bakelite resin for ease of handling, a section was cut 
perpendicular to the plane of fracture. 


* Crown copyright reserved. Reproduced with the permission of the Controller of H.M. 
Stationery Office. 

+ Slip-casting consists of pouring finely grcund solid particles suspended in a suitable liquid 
medium into a porous mculd, in which a cast is formed when the liquid is absorbed. 
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Plate 1. 


Plate II. 


Plate III. 


Fracture in a sintered aggregate of corundum crystals. (Magnification: 100x 3) 
To face page 249 
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‘The section was ground with silicon carbide particles of successively finer 
grades on a wet glass plate. ‘The finest grit had a maximum particle size of 12 p. 
‘The section was then polished in two stages; firstly, boron carbide (maximum 
particle size 6 ..) was employed on a pitch lap and, secondly, diamond dust (particle 
size range 4 to <1) on a lead lap. 

The polished section was etched by immersion in molten potassium hydrogen 
fluoride for approximately one minute. 

The above procedure allowed study of that part of the fracture intersecting a 
plane cut through the interior of the bar. It was also desired to know the mode of 
fracture at the surface, for which purpose some bars were polished on their 
tension and compression faces, fractured, and etched afterwards. 

The prepared sections were examined with the aid of an orthodox reflected 
light microscope using unpolarized light. The photomicrographs in Plates 
I-III are 100 times natural size. 

The bars were made by slip-casting, the particle size being <10y prior to 
sintering at 1,650°c. for 4 hours. The specific gravity was 3-80. 


$3. RESULTS AND DISCUSSION 


Plates I-III show interior sections of fractures. ‘Transcrystalline failure has 
occurred at A-A, B-B, E-E. Shattering has caused removal of much material 
from other places adjacent to the fracture. The angular outlines of C, D (Plate IT) 
suggest that they are parts of broken crystals. It is also considered that inter- 
crystalline failure occurs, but evidence for this type is necessarily less definite than 
for transcrystalline failure. Study of twenty fractures strongly suggests that 
intercrystalline failure occurs in some cases. For instance, the irregular shape of 
the exposed perimeter of crystal F (Plate III) is more like a crystal-boundary than 
the types of transcrystalline failure observed. 

No difference has been detected between the types of fracture occurring at 
the surface of and inside a bar. 

Two types of transcrystalline failure may be distinguished; some fracture 
surfaces appear to be accurately planar, while others are curved. This is con- 
sidered to depend upon whether the crystal has failed along a parting plane or has 
- developed conchoidal fracture. Most writers on corundum agree that there is no 
true cleavage, but that crystals may sometimes be broken readily along {0001} 
basal and {1011} rhombohedral planes (corundum is trigonal, class 3m). The 
phenomenon is called “‘ parting”; in the case of the rhombohedral planes it is 
often attributed to polysynthetic twinning. If parting does not take place, 
corundum is said to fracture conchoidally. Crystals A-A, B-B, C, D appear to 
have failed conchoidally, whilst E-E appears to have parted. 

Crystal-boundaries do not all show at the same stage of etching, so that when all 
are uncovered there is a divergence in the extents of attack at the different bound- 
aries, as can be seen in the photomicrographs. Furthermore, the extent of attack 
frequently varies along the length of curved boundaries. It is considered that the 
extent to which a given boundary is attacked is a function of the difference between 
the orientations of the crystals forming the boundary and a function of the direction 
of the boundary. There is a connection here with recent work in metallurgy by 
Lacombe and Beaujard (1947), Lacombe (1948), and Forsyth, Metcalfe, King and 
Chalmers (1946), who have shown that the relative orientations of a grain boundary 
and the adjacent grains influence a number of phenomena, including etching, 


250 | ¥. P. Roberts 


boundary precipitation and ‘thermal etching”. These studies are closely 
related to the grain-boundary transition-structure hypothesis of Gough (1927), 
Hargreaves and Hills (1929) and Chalmers (1937). 

Well-formed etch-figures may be produced in the individual corundum 
crystals, provided that a planar polish of high quality is obtained and extreme care 
is taken during etching. A detailed study of the figures would no doubt provide 
useful information on the crystallographic orientations of the crystals which, in 
turn, would probably allow a closer investigation of the above-mentioned differ- 
ential etching and perhaps lead to information about the mechanical strength of 
different types of boundary. 

It is thought that the occurrence of trans- and intercrystalline failure at each 
pcint in a fracture face has been decided by one or both of two factors. Corundum 
crystals are anisotropic. The weaker planes within the crystal probably undergo 
fracture easily, whilst the stronger planes may be sufficiently strong to deflect the 
course of a fracture round the crystal-boundary. In this way the anisotropy of the 
crystals could account for the occurrence of simultaneous trans- and inter- 
crystalline failure. A second factor which also would account for the experimental 
results is the non-uniformity of the boundaries, to which reference has already been 
made. 

It is believed that reflected light examination of sintered corundum and other 
similar materials will prove to be superior to transmitted light examination for 
many studies, particularly those associated with crystal-boundaries. The author 
has used the technique to advantage in a study of bodies with high porosity. 


§4. CONCLUSIONS 

Since both trans- and intercrystalline failures occur when bars of well-sintered 
corundum are broken, it is concluded that the strengths of the crystals and the 
crystal-boundaries must be of comparable magnitudes. The results of experiments. 
performed on large corundum single crystals, made by the Verneuil process and 
machined into bars, do not conflict with this conclusion. ‘The bend strength - 
of a single crystal is approximately equal to that of a well-sintered polycrystalline 
similar bar, suggesting that the crystal-boundaries are either as strong as or stronger 
than the crystals. 

Reflected light examination of sintered corundum is not only possible but is in 
many respects superior to transmitted light examination, as, for example, in this 
type of investigation. 
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ABSTRACT. A torsion device is described in which a small and decreasing stress is 
exerted on the visco-plastic substance by means of a blade attached to a torsion wire. A 
modification of the Scott Blair-Nutting equation is suggested to account for the variation of 
both stress and strain with time and a “‘ firmness intensity factor’ of the substance is 
calculated, by means of the modified equation, from the observed relation between blade 
deflection and time. The apparatus has been used to study the properties of milk curds. 


Si ENC ROD U CirowN 

UCH work has been published on the elasticity of gels, but little informa- 
M tion is available on the region of plastic flow. The latter is particularly 

important for weak gels, where the elastic limit is easily exceeded, and 
where, therefore, the strength of the gel is primarily dependent on its visco- 
plastic properties. Measurements of such properties have been carried out 
hitherto on an empirical basis by cutting the gel by means of a specified knife 
or a plunger (Bloom 1924, Miller 1935) and thus subjecting it to an intensive 
local stress. This method has obvious limitations. Measurements under 
conditions of constant applied stress, which have been made on a number of 
plastic materials (Scott Blair and Coppen 1941, Broome and Bilmes 1941), and 
measurements of the shear force required to fracture gels (Hamer 1947), were 
found to present considerable experimental difficulties in the case of weak gels. 
Cheftel and Mocquard (1947) developed an apparatus for the recording of the 
deformation of gels under a varying applied force, but it did not permit evaluation 
of specific rheological properties. As a preliminary to the investigation of the 
visco-plastic properties of weak gels, therefore, a method of measuring their 
rheological properties was developed which is simple and free from some of the 
theoretical objections of the arbitrary methods. ‘The method also permits 
determination of the viscosity of Newtonian liquids, and of the elasticity of 
Hookian solids, which may be considered as limiting cases of the visco-plastic 
state. 

5 2 AP PAR AMINES 

The rigidity of a material is found by observing the displacement under an 
applied force. ‘The force should be exerted over an appreciable area, so that 
variations in the homogeneity of the materials are averaged out, and excessive 
local stresses are avoided. The most convenient way of applying a small stress 
is by means of a comparatively large blade attached to a fine torsion wire. 

The construction of such a torsion apparatus is shown in Figure 1. A brass 
blade B is attached to a torsion wire W by means of a suitable clamping device. 
The wire is held by another clamp fixed to the torsion head H. ‘Two pointers U 
and L are attached to the upper and lower clamps respectively. ‘The apparatus 
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is mounted on a retort stand so that the blade B can dip into a dish containing the 
substance, the firmness of which is to be measured. ‘This dish is placed in a 
thermostatically controlled water bath. 

If the upper pointer U is moved rapidly through a given angle between two 
adjustable stops, by turning the torsion head H, the blade B turns in the substance, 
and its deflection can be found by noting the position of pointer L on the scale S 
at given intervals of time. The rheological constants of the material can then be 
calculated from the relation between pointer deflection (in degrees) and time 
(in seconds), as described below. 


Maximum Deflection 


ewe —_ —— 


b 
= 


Deflection & (degrees) 


0 50 100 150 200 
Time f (seconds) 


Figure 1. Diagram of torsion device. Figure 2. Relation between blade deflection (8) and 
Limiliorsionihead: time (t) for oil and milk curds. 
U. Upper pointer. 1. Polymerized castor oil (y=120 poises). 
W. Torsion wire. 2. Boiled milk (pH 5-02). 
L. Lower pointer. 3. Full cream milk food (pH 5:3). 
B. Blade. 4. Pasteurized milk (pH 0°52). 
S. Scale. 5. Full cream milk food (pH 5:86). 
D. Glass dish. 6. Dried skimmed milk (pH 5-95). 


The dimensions of the blade used were 1}in. x 3in.x tin. The blade 
thickness is immaterial as long as the weight is sufficient to keep the torsion wire 
taut, but the inertia of the blade must be small compared with the retarding force 
exerted by the medium. If the diameter of the dish is not less than twice the 
length of the blade, no slipping of the substance occurs when the blade is turning. 
"Tungsten and beryllium—copper are suitable metals for the torsion wire. The 
length and diameter of the wire depend on the maximum stress required. A 
6in. length of 40s.w.c. tungsten wire, giving a maximum stress of 72 dynes/cm? 
with the 14in. blade for torsion head deflections of 60°, was found to be suitable 
for the measurement of viscosities of the order of 100 poises. 


§3. EXPERIMENTAL WORK 
The present application of the torsion device was concerned primarily with 
the investigation of the properties of milk curd. This curd forms a weak gel 
with a low elastic limit. ‘The properties of the gel vary with the acidity and the 
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‘type of the milk, and a wide range of different curds could thus be obtained. The 
curds were formed by coagulating an acidified 12:5°% solution of milk powder 
by means of 0-45°% solution of pepsin at 35°c. The torsion blade was inserted 
before the formation of the curd. Cutting by the blade destroys the structure 
of the curd, and a test cannot, therefore, be repeated on the same sample. This 
imposes some limitations on the reproducibility of the deflection-time curves 
where random variations from sample to sample occur, as was found to be the 
case with milk curds. 

Curves of blade deflection against time for a vegetable oil and different milk 
curds are given in Figure 2, the maximum (torsion head) deflection being 60°. 
The readings of deflection were taken at intervals of ten seconds up to two minutes. 
Accurate readings could not be obtained below ten seconds because of the rapid 
initial movement of the blade. Some errors are also to be expected at the 
beginning of the deflection-time curve, owing to the impact of the blade on the 
medium. 

§4. THEORETICAL TREATMENT 

The deflection-time, 5-t, curve for oil, which is a Newtonian liquid, 
approaches the maximum deflection, i.e. the torsion head deflection «, logar- 
ithmically, and it can thus be expressed by the exponential relation 

Cem OXD A Ab) oy Ch ene LE (1) 
where A is a constant. 

The curves for milk curd show that the deformation of these materials can 
be expressed neither independently of time, as in the case of Hookian solids, 
nor as an exponential function of time, which would indicate the viscous flow 
of a Newtonian liquid, and we are, therefore, dealing with ‘“ intermediate” 
materials of negligibly small elasticity. 

Nutting (1921) and Scott Blair and Coppen (1940) have suggested an equation 
relating the shear strain o of such materials to a power f of a constant shear stress S 
and also to a power & of the time of application ¢ of the constant stress, so that 

Ga aie, te ee ate (2) 
where ’ is a constant of the material which may be described as the “ firmness 
_ intensity factor”’.* This equation is not valid if both o and S vary with ¢. 

For the torsion device described above, S is not constant but varies with o 
 andthuswith¢. Several equations can be established to account for the variation 
of c and S with ¢. Because of its simplicity, Scott Blair’s equation for the rate 
of straining (Scott Blair 1944) is suggested in which the rate of change of o with 
tis given by 

AGE ee cyto 8 ene (3) 
Sj is the mean stress over the short time dt and is assumed constant during this 
time. This equation is applicable to the torsion device only, where the variation 
of o and S with ¢ are interdependent. It is equivalent to equation (2) when 
Sq 1S constant. 

For a viscous liquid, B=1 and k=1 and hence 
ANCIENTS OS Nel aagilel Bae Seo (4) 

Since doidt aay dx) dt=avjadx® rr i Cae elles (5) 
the velocity gradient dv/dx can be substituted in equation (4), which then becomes 


* This term is used at the suggestion of Dr. G. W. Scott Blair. 
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Newton’s equation of viscous flow with ‘Y equal to the viscosity 7. For an 
elastic solid, B=1 and k=0, so that do/dt=0. This represents the equilibrium 


condition between applied stress and reaction which will be set up almost instan- 
taneously by the torsion blade, and is independent of time. Equation (3) is 


~~ 


thus shown to be valid for two particular cases, and it may, therefore, be taken to _ 


be valid for any case intermediate between the elastic (solid) and the viscous 
(liquid) state. 

Equation (3) can only be solved generally for o if Sw is specified. For a 
torsion wire of diameter D and length L twisted through an angle @ radians, the 
torque T is given by 

TenODEn 321s wre qintanh Tl eee (6) 
where N is the rigidity modulus of the material of which the wire is made. 

The torque per degree twist, which is a constant for a particular wire, is there- 
fore 


G=anT/1800=n2D*N/S760L. nae (7) 


If the torsion head is turned through an angle «, the torque at a blade deflection 
x 1S 


T= GG 8) et ee eee (8) 
and the corresponding stress exerted on the substance to be tested is 
Sy=4G(«—a)la- = 2" Sheth Bee (9) 


where / is the length of the blade and a the total blade area. For a particular 
wire and blade arrangement, 
Sip = Cy(x — x) 
C, being a constant for the instrument. 
It is more difficult to define the strain as there is certain to be some flow at 
the edges of the blade as the blade moves. But it is reasonable to define an element 
of strain as being proportional to an element of deflection, and hence 


do.=Gedx: Pte Ee (11) 
Substituting (10) and (11) in equation (3) and replacing the parameter ¢ by the. 
variable ¢,, we obtain 
dx|(a —x)P = Cekt*—1dt, [¥C,. 
Integrating between the limits x=0, x=6, t,=0, t,=t, where 3 is the blade: 
deflection at time ¢, this becomes 


(0) So O(a Cr TC, -2 1 eee (13) 
Equation (13) can be simplified when 5/« <1, so that 
devdla\ Vet (Bal) Sieve ne apr eesti (14) 


and, therefore, from equation (13), 


S= (aC, EC: 


Writing, «C;=C, where C is an instrument constant, the “firmness intensity 


factor” ‘V’ is given by 
P= CereG.. 
Equation (16) is expressed in logarithmic form as 


logd Slog (C8PC,) Hlont pay anes (17) 


Rheological Properties of Visco-plastic Substances 255, 


__ The relation between logd and logt should, therefore be linear. k, which 
is usually called the ‘dissipation coefficient’’, is given by the slope of this line. 
B is found for a given material by plotting log 5 —logt lines for two experimental 
conditions, C’ and C”. If the corresponding intercepts with the logé axis 
are B’ and B’, then 


Dena ioe (CHG a Se (18) 
If 8 =1, integration of equation (12) results in an equation 
Inga (@— oN CeCe Hts (19) 
or, taking logarithms and substituting Cy, =C,/C,, 
loginia/(a—o)i=log(Ce/ Per logit. 2 =- a... (20) 


The linear relation between log In («/(«—8)) and log¢ for the vegetable oil and. 
for some milk curds is shown in Figure 3. The curves become inaccurate at. 
short times for the reasons discussed above. 
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Figure 3. Relation between log In a/(a—8) and log ¢ for oil and milk curds. 
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. Polymerized castor oil. 

. Full cream milk food (pu 5:02). 
. Boiled milk (pH 5-02). 

. Modified milk food (px 5-07). 

. Raw milk (pH 5:53). 

. Pasteurized milk (pH 6:0). 


If 8/a<1, In {x/(« —8)} =5/a, and equation (16) can also be applied to any 
material for which B&1. The use of a separate equation for the special case 
when B=1 can thus be avoided. The relation between logs and logt can be 
expected to be linear where the relation between log In(«/(«—6)) and logt is 
linear. The linearity of the curves given in Figure 3 is not dependent on the 
value of f. 

For viscous liquids, k=1 and ‘ becomes the viscosity 7. Equation (19), 
can then be simplified, so that 


hia: ) SO (21) 
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If this equation is expressed as an exponential relation, it becomes 
E—~O=a expt (Cre Ne eee (22) 

This is the same as equation (1), with A=Cy/n. It is thus shown that equation (3) 
is the general equation applicable where both o and S vary with ¢, as is the case 
‘in the torsion apparatus described above. 

For an elastic solid, k=0 and VY becomes the modulus of elasticity E. 

Hence, 

lnfol(d@—Sih= CyB 9 5 9 andtes (23) 

Thus, the “ firmness intensity factor’ ‘ of a visco-plastic material can be 
found, where 81, from equation (16), and, where B=1, from equation (19). 
It should be emphasized that in order to compare different materials it is not 
sufficient to give only values of ‘’. The values of k and 8 depend on the 
material and should be given in addition. The value of k, which lies between 0 and 
1, indicates how far the material departs from a true liquid or from an elastic solid. 
Equation (16) can also be used where B=1 if 5/a«1. The viscosity 7 of a 
Newtonian liquid is given by equation (22). If a material has an appreciable 
‘elasticity, the modulus of elasticity E can be computed, by means of equation (23), 
from the deflection at the equilibrium position of the torsion blade for a given 
applied stress. 

§5. CALIBRATION 


The torsion device is calibrated by means of a liquid of known viscosity 7. 
If a log In {a/(a—8)}-log ¢ curve is plotted for such a liquid, the ratio Cy/7 is 
given by the intercept of the line with the log In {x/(« —5)} axis. Cy can then be 
-calculated and used to find unknown values of ‘ or 7 from equations (19) or (22). 

C, has to be determined before ‘’ can be calculated from equation (16). 
‘This can be done by assuming that C, is not affected by changes of the constants 
‘, B and n, so that C, can be derived from tests on a different medium for which 
‘these constants are known. Since C, can be calculated, C, can be found from Cy. - 
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ABSTRACT. Four methods of conditioning a sample of material to a prescribed liquid 
or vapour content are examined mathematically. The ‘“‘ two-stage method ” of conditioning 
leads to considerably shorter conditioning times than the other methods provided it is carried 
out in a certain way, and information presented in graphical form in this paper shows how 
to obtain the maximum advantage in practice. The necessary equations are developed for 
both a plane sheet and a cylinder, and the corresponding problem in heat flow is mentioned. 


SN ER ODT Ce iiON 


COMMON requirement, both in laboratory and industrial practice, is to- 

cause a sample of material to absorb a prescribed amount of a given liquid 

or vapour, a process usually referred to as the “conditioning” of the 
sample. ‘This frequently means that the sample is conditioned to be practically 
in equilibrium as regards its total content with a vapour of a given pressure, but 
there is sometimes the additional requirement that the distribution of liquid or 
vapour through the sample must have reached a prescribed degree of uniformity. 
The corresponding problem in heat flow is concerned with the conditioning of the 
sample to a uniform temperature. ‘The problems are closely allied in that each is 
essentially a problem in diffusion, the first diffusion of matter and the second 
diffusion of heat, and the basic mathematical equations are the same in the two 
cases. ‘The most obvious and simplest method of conditioning a sample to a 
vapour of given pressure is to suspend the sample in a bath or enclosure containing 
the vapour until the sample reaches equilibrium, the vapour pressure being 
maintained constant at the required value. In some cases, however, for example 
the conditioning of polymer films to water vapour or organic vapours at low 
pressures, the process is a slow one, and it is not uncommon for the time to condition 
satisfactorily to be as much as three or four weeks. It 1s therefore worth while to 
explore possible ways of reducing the time of conditioning. ‘The present paper 
is concerned with a theoretical study of four ways of carrying out the conditioning 
process, and in particular with the calculation of optimum conditioning times for an 
ideal case in which conditioning is considered as a simple diffusion process in which 
the diffusion coefficient (or heat diffusivity) is assumed constant. Later in the 
paper, one example, in which the diffusion coefficient is assumed to increase 
linearly with concentration of diffusing substance, is included to indicate the 
general effect of non-constant diffusion or thermal properties and the extent to 
which the ideal calculations may be used as a guide to experimental procedure in a 
practical case. 

The present discussion relates to diffusion of matter because the incentive to 
carry out the calculations described came from experimental work on diffusion in 
polymer films. The equations and general conclusions apply equally well to heat 
flow, though in that case uniform temperature distribution is nearly always the 
criterion of conditioning rather than the attaining of a prescribed total heat content. 
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§2. ALTERNATIVE METHODS OF CONDITIONING 
Four different methods of carrying out the conditioning process are examined in 
this paper. ‘They may be classified conveniently as (a) infinite bath conditioning ; 
(b) two-stage conditioning ; (c) finite bath conditioning ; (d) insulated conditioning. 
‘The significance of each term needs some explanation. 


(a) Infinite Bath Conditioning 


This is the straightforward method referred to above in which the sample to be 
conditioned is immersed in a bath of vapour, the vapour pressure in the bath 
remaining constant at the value to which the sample is to be conditioned. ‘The 
condition of constant vapour pressure is most readily secured by having a large 
amount of vapour in the bath compared with the amount finally taken up by the 
sample; theoretically this implies an infinite bath of vapour. 


(b) Two-stage Conditioning 

In this method two infinite baths of vapour are used, one at the pressure to 
which the sample is to be conditioned and the other at a higher vapour pressure. 
For convenience in discussion this will be referred to as saturated vapour though 
in practice it will be as high as practical considerations allow. ‘The sample is first 
placed in the saturated vapour and after a certain time is transferred to the other 
bath and left to attain equilibrium. It will be shown that the success of this 
method depends on the time at which the transfer is carried out. 


(c) Finite Bath Conditioning 


If there is only a restricted amount of vapour, the vapour pressure in the bath 
decreases steadily as vapour is absorbed by the sample, the rate of fall being 
determined by the condition that the total amount of vapour in the bath and in the 
sample is conserved. If, for example, the vapour in the bath is initially saturated, 
the pressure can be made to approach the value to which the sample is to be 
conditioned by suitable adjustment of the respective volumes of bath and sample. _ 


(d) Insulated Conditioning 


The sample is again immersed in an infinite bath containing vapour at the 
highest convenient vapour pressure, and when the total content of the sample 
reaches the prescribed value the sample is removed from the bath and insulated by 
some means so that no further change in total content takes place while the 
desired degree of uniformity of concentration through the sample is approached. 
It will be seen that the same result can be achieved by causing the vapour pressure 
in the bath to vary in a particular way instead of insulating the sample. 


§3. MATHEMATICAL SOLUTIONS FOR A PLANE SHEET 

Subject to certain simplifying assumptions, calculations of the rate of condition- 
ing by each of the above methods can be made. In this and the following two 
sections the sample is considered to have a geometric form such that it approxi- 
mates to a solid of finite thickness bounded by two planes, e.g. a disc the area of 
whose edge is negligible compared with the area of the plane faces, so that condition- 
ing may be considered to take place by diffusion in one dimension only, normal to 
the plane faces. ‘The diffusion coefficient (or heat diffusivity) is assumed constant. 
In treating the finite bath problem it is also necessary to assume Henry’s law of 
solubility in order to obtain formal solutions of the mathematical equations. For 
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} conditioning by the other methods, a knowledge of the absorption isotherm may be 


} necessary in order to apply the methods in a practical case but it is not necessary for 
| the mathematical treatment presented here, and the isotherm need not be a linear 
| one. Since the liquid or vapour to which the sample is to be conditioned may be 
| either a complete solvent or only a swelling agent for the sample, it is convenient 
to refer to the diffusing substance by the general term “ penetrant”’. 
If C is the concentration of penetrant within the sample and D the constant 
} diffusion coefficient, the basic equation, of which solutions are required for the 
boundary conditions corresponding to each of the four methods of conditioning, is 
2 
eo 8) ae See $ Zeal a dia oils pid (1) 
ot thee 


In the corresponding equation for heat flow C is replaced by temperature and D 
becomes the heat diffusivity given by x/ps where « is heat conductivity, p the density, 
and s the specific heat of the sample. When the sample is immersed in an infinite 
bath it is supposed that the two surfaces of the sample reach a state of equilibrium 
with the vapour in the bath instantaneously, the concentration of penetrant just 
within each surface being related to the vapour pressure in the bath by the absorp- 
tion isotherm. If the surface concentration of penetrant is denoted by Co, a 
solution of (1) is required for the boundary conditions 


C=C; Vel, x=, Joy o Bien: (4) 
and the initial condition 
C=0; Ozen <al, BeOS warn \Pe eet Tos (3) 


J being the thickness of the laminar sample. The solution is well known (Carslaw 
and Jaeger 1947) and is 


1. (2m+1)rx { — D(2m+ | 
n-——_—— exp , ———_,, 


Sl RP 


2 | Y 4 = 
€=0,|1-5 3 5 l 


The total amount of penetrant M, absorbed by the sample at time ¢ is given by 


l 
integrating (4) with respect to x, i.e. M, =| C(x’, t) dx’, and it is easily shown that 
0 


8 1 — D(2m + 1)?x?t) 
M=Cd| 1= G3 ge apexP{—— Ft I Paes (5) 
The first part of two-stage conditioning, prior to transfer to the second bath, is also 
governed by equations (4) and (5), with C, substituted for Cy if C, is the constant 
concentration of penetrant just within the surfaces of the sample in the first bath. 
At the time of transfer, say t,, there will be some sort of non-uniform distribution of 
concentration through the sample so that diffusion in the second bath is subject to 
conditions of the form 


C=C, =constant, A). Meat, Fie! We da (6) 
Cai x), One xt c= SID Jar, (7) 
and in general the solution of (1) for conditions (6), (7), is 


= Qa 2( p 
1 sin Cet D7 0 { Acct as all] 
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This can be simplified in the present context because f(x) represents the concen- 
tration distribution in the sample at the end of the first stage of two-stage condition- 
ing, i.e. it is the result of immersion in an infinite bath of concentration C, for a time 
t,. f(x) is thus given by substituting the values t=¢,, Cy=Cyin (4). In this case 
it is easy to show that 
y . Mas ZiCs — D(2r +1)?x?t, ] 
’) sin —— dx’ = ~——~| 1- <> |, ove 9 
| fe ) sin 7 dx (r+) exp { p } (9) 
where (27 + 1) is written for odd 7, so that the course of conditioning in the second. 
stage is determined e the equation 
. (2r+-1)ax — D(2r + 1)?x(t — a} ] 
C=C, ce me E Sir ike) pete pel ges ale 
AG aed ‘1 (or +1)rx 
TT r=0 2r + i i 


a fe eae cal ne RS cy 


[2 
The integrated form of (10) is 


$ ea utes) scl 21 
M=Cal[1— © Gao P +t 2, rte 
an 22 rary 

x [ 1 exp {POE exp {ee wee (11) 


The problem of the finite bath may be formulated mathematically as follows. 
Suppose the sample, of thickness /, occupies the space — 3/<x<$l, while the bath, 
of limited extent, occupies the spaces —3I +h) <x<— 41, sl<x<h(l+h). The 
concentration in the sample is initially zero, and in the bath is always uniform. _ 
Let the partition factor k denote the ratio of the concentration in the sample to that 
in the bath when a state of equilibrium is reached, and C, denote the initial surface 
concentration just within the sample. ‘l’hen the initial uniform concentration in 
the bath is C,/k. It is assumed that the surface concentration just within the 
sample is always k times the uniform concentration in the bath throughout the 
conditioning. ‘Taking C(x, t) for the concentration of penetrant in the sample as 


“41 
before, the total amount of penetrant in the sample at time ¢ is 2 C(xt,.2) dx. 
0 


and so the boundary condition on the surface x =4/, expressing conservation of 
penetrant is 


41 
hC|k+2 | CdeanCik . aid atten meaeeeeen (12) 
0 


The initial condition is 
Ca —H<x<chl, tai fal payee rk. (13) 


A convenient solution of the diffusion equation, subject to conditions (12), (13); 
takes the form . 

© ’ 1+a  cos(2q,x/) 

= I ae) 3; a Dg?) ae ee 

C a + eRe D4;,t) [power an casa: }, owiehe (14) 
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| where the q,,’s are the non-zero positive roots of tang, = —«q,,; «, the parameter 
expressing the finiteness of the bath, is given by «=h/(kl); C,, is the uniform 
concentration in the sample attained theoretically after infinite time and 


C,,=4C,/(1+«). The integrated form of (14) is 


M,_,_ & 2x(a+1) 
M., =] Penseeeee et ao) sioteterere (15) 


where M,, is the amount of Lean ar, in the sample after infinite time and is 


related to ie other parameters by M,,=C,h/k(1+«). The solution (15) was 


derived by A. H. Wilson (1948) and (14) follows readily from the equations of his. 
aper. 

; The process of conditioning by the fourth method prior to insulation is again. 

governed by equations (4) and (5) for the infinite bath. The condition that the 

sample is insulated after a certain time ¢, is simply that dC/dx =0, x =0, «=I, tty. 

The solution of the diffusion equation for a plane sheet insulated on both surfaces 

and with the initial condition C =f(x), 0 <x </, t=t,, is (Carslaw and Jaeger 1947) 


rl (oe) AF Biel Fas 
wil ey peu DEX srevandba ts) | 03 | He) Led pp i 
LJo dja ie l 


In the present example f(x) is given by evaluation of (4) at the time t=¢,, with C, 
written for C) if the infinite bath used prior to insulation is taken to have the same 
concentration as that used in the first stage of two-stage conditioning. It can 
readily be shown that the two integrals of (16) then reduce respectively to 


= 1 — D(2m + 1)?7t, 
i; fla’) dx’ =c1/1-5 2 Gmaiye xp eh, pa (17) 
and 
(ae 2prx' ,, _ 8 les 1 —D(2in-P 1776, | 
| fee ) cos jj dx = 2 ; 4p? — (2m + 1 exp{ 5) 5 
hopdere (18) 


where 29 is substituted for 7 since only the terms involving even m are non-zero. 
The property of particular interest is the variation of surface concentration with 
time after insulation, and this is given by evaluation of (16) at x=0, having sub- 
stituted from (17) and (18) for the integrals. We have finally, 


Pena Sl gp f= an Uh) 
Gili a ci a (2m +1)? es ip 
162 2 1 silat { — D(2m +1)2n°t, 
4 ee Se " J bese [2 ° 


Equation (19) expresses the way in which the concentration at the surface of the 
sheet varies after insulation, the total content of the sample remaining unchanged. 
As an alternative to insulation if the surface concentration is caused to vary by 
some external means, in accordance with (19), after the time ¢,, the effect will be 
identical with insulation, the surface gradient of concentration remaining constant 
and zero, and the sample neither gaining nor losing penetrant. 
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§4. NUMERICAL COMPARISON OF THE FOUR METHODS 
OF CONDITIONING 


In order to see whether the time of conditioning differs appreciably in the 
different methods, some of the above solutions have been evaluated numerically. 
Since it is only after infinite time that a sample is completely conditioned to the 
extent of having a perfectly uniform concentration throughout, some arbitrary _ 
criterion of a conditioned sample must be adopted. A sample will be said to be 
conditioned with respect to its total content at a time when the content differs by 
+1°%, from the ideal value it would reach after infinite time, provided it never 
deviates by more than + 1° from the ideal value at any subsequent time. The 
distribution of concentration through the sample at the time of conditioning, as just 
defined, is examined later. 

The constant concentration just within the surface of the sheet in the first part 
of two-stage conditioning has been denoted by C, and it is convenient to use this as a 
unit of concentration. The corresponding total content of the sample, after 
immersion for an infinite time in an infinite bath such that the surface concentration — 
is C,, is 1C,, which is a convenient unit in which to measure total content of pene- 
trant inthe sample at anytime. Thus on this scale, the total content of the sample 
will approach unity if left in the first stage bath of two-stage conditioning for an 
infinite time. The non-dimensional unit of time Dt//? is used as usual. Figure 1 
shows the progress of conditioning to a total content of 4C, in an infinite 
bath for which the surface concentration is 4C,, and also by two-stage 
conditioning. ‘The curve shown for the infinite bath is obtained by evaluating (5) 
for Cy =4C,, and the curves for two-stage conditioning by evaluating (5) for Cy) =C, 
and (11) for C,=4C,. A number of curves are drawn corresponding to different 
times of transfer ¢, from the first bath to the second. ‘The four values of time of 
transfer are given by Dt,/]?=0-012, 0-025, 0-040, 0-060 respectively. The curve 
corresponding to transfer after the shortest time shows a minimum, the total 
content of the sample first decreasing after transfer and then increasing and 
approaching equilibrium at a time not much different from that shown by the 
infinite bath solution. The curve for which D¢,/l? =0-025 has a similar behaviour. 
‘The sample left in the first bath for the longest time contains an excess of penetrant 
at the time of transfer, and after transfer the total content steadily decreases reaching 
the 101% limit only a short time before the infinite bath solution reaches the 99% 
limit. Figure 1 shows, however, that there is an optimum time of transfer at about 
Dt,/l? =0-040, for which the conditioning curve reaches the 101°% limit much 
sooner than any of the others, and which remains within the limits of +1% 
thereafter. This optimum curve has a minimum at 99°% of the equilibrium 
content. In Figure 2 the optimum curve for two-stage conditioning is compared 
with the infinite bath curve and with the finite bath curve calculated from equation 
(15) with « = 4, the initial surface concentration in the finite bath being C,, the same 
as in the first bath for two-stage conditioning. The finite bath method leads to 
quicker conditioning than the infinite bath, but is slower than the optimum two- 
stage conditioning. ‘Thus the times of conditioning by two stage, finite bath and 
infinite bath methods in this example are given by Dt/l?=0-062, 0-200, 0-450 
respectively. If the two-stage conditioning is prolonged till the 99° limit is 
reached, to be more strictly comparable with the other two, the time 0-062 is to be 
replaced by 0-090, so that two-stage conditioning represents an improvement on 
infinite bath conditioning by a factor of five or more in this particular case. 
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another at a different concentration in two-stage conditioning would produce a 
non-uniform distribution of concentration within the sample which may persist 
from some time after the transfer, so that, even when the sample is satisfactorily 
conditioned as regards its total content, it may not have achieved the desired 
degree of uniformity of concentration. Figure 3 shows concentration~distance 
curves obtained by evaluating (10) for C,=4C, for the times at which M/C, 


: It might be expected that the transferring of the sample from one bath to 
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Figure 3. Concentration—distance 
curves for a plane sheet. 
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Figure 4. Conditioning and trans- 
fer times for a plane sheet. 
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Figure 2. Comparison of methods of conditioning of B. Time to reach 99% equilibrium 
a plane sheet Cp=}Cj. content. 


reaches 101% and 99% of the final equilibrium value of $/C,. At the earlier time 
the concentration is appreciably non-uniform, the maximum deviation from the 
desired concentration of C/C,=4 being about 10%. By the time the total content 
has fallen to the 99° limit, however, the concentration is as uniformly distributed 
as it would be after conditioning to 99°% in an infinite bath, the concentration 
distance curves being indistinguishable on the diagram. 
; 18=2 
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Corresponding calculations for conditioning to different total contents lead to 


qualitatively similar results as regards the overall rate of absorption curves. 
In each case there is one time of transfer in the two-stage method which leads to a 
substantially shorter time of conditioning than is achieved by finite or infinite bath 


; 


| 


methods. It is not possible to show all the results in detail but in Figure 3 


concentration—distance curves are shown for cases in which the total content of the 
sample is conditioned to 2/C,/3 and 9/C,/10 in addition to 3/C, already discussed. 
When the total content reaches 101° of 2/C,/3, the concentration shows a maxi- 
mum deviation of about 14° from the desired uniform concentration of 2C,/3, 
but when the 99°, limit is reached the corresponding maximum deviation of 
concentration is only about 3%. The concentration—distance curve is also shown 
for the time at which the total content reaches 99°, of 2/C,/3 by finite bath con- 
ditioning. For the case in which the equilibrium total content is 9/C,/10, the 
101% limit is never attained in two-stage conditioning, the value of M,//C, at which 
the sample is transferred to the second bath lying within the +1% limits. In this 
case when the 99% limit is reached the concentration distribution is seen to be as. 
uniform as possible. Inspection of concentration—distance curves, evaluated for 
intermediate values of equilibrium total content but not shown here, reveals that 


the maximum deviations of the concentration from uniformity when the 101% and © 


99° limits are reached are most serious for the case in which the equilibrium total _ 


content is 2/C,/3; the deviations decrease for higher and lower total contents and 


so in two-stage conditioning, when the sample is conditioned to 99°% of the © 
equilibrium total content, the greatest deviation of concentration at any point in the — 


sample from the desired uniform value is never greater than about 3%. 


In Figure 4, the lower curve shows the way in which the optimum time of 


transfer from the first bath to the second in two-stage conditioning depends on the 
equilibrium total content desired, i.e. the concentration in the first bath is always 
taken to be C, and in the second bath the concentration takes different values 
depending on the prescribed state of the conditioned sample. ‘The upper curve 
shows the corresponding behaviour of the time taken to condition to 99% of the 
equilibrium total content. Both times increase as the equilibrium total content is. 
increased. In practice this means that a given equilibrium content will be 
achieved most rapidly by making the concentration in the first bath as high as 
possible subject to whatever practical limitations may exist. ‘This is obvious. 
qualitatively since the higher the concentration in the first stage the greater the 
concentration gradients involved and hence the quicker the diffusion. To 
estimate the effect quantitatively, however, necessitates the mathematical analysis. 
leading to Figure 4. For a given equilibrium content the time of transfer may be 
read from the lower curve provided the diffusion coefficient D and thickness of the 
sample are known or can be found for the system concerned. Both curves reach 
the value Dt//?=0-45 when the equilibrium content is /C, since we then have 
conditioning by the infinite bath method arising as a special case of two-stage 
conditioning. ‘The effect of making the concentration in the first bath high 
relative to that in the second bath is also to produce a more uniform distribution of 
concentration through the sample at the time of conditioning as previously shown 
by Figure 3, provided the desired equilibrium content is less than 2/C,/3. 

For practical purposes in two-stage conditioning it may sometimes be more 
convenient to know the value of the total content of the sample at which the transfer 
is to be made rather than the time of transfer, i.e. to know M,//C, at time ¢, instead 
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of Z, itself. In Figure 5 the excess absorption at the time of transfer is plotted as a 
; function of the equilibrium total content. Excess absorption is the amount by 
_ which the total content of the sample at the time of transfer exceeds the desired 
) final equilibrium content. This excess is necessarily zero at each end of the curve, 
| i.e. for equilibrium content zero and /C,, and the curve passes through a maximum 
_ value for an equilibrium content of about 0:16/C,, where the excess absorption is 
0-146/C,. For equilibrium contents higher than that at which the maximum 
occurs the relation plotted in Figure 5 is closely linear. Thus if the total amount of 
_ penetrant in the sample can be measured by some means as conditioning proceeds, 
_ Figure 5 may be used to determine, for a prescribed equilibrium content, what the 
_ content of the sample should be when it is transferred from the first to the second 
bath in order to secure optimum two-stage conditioning. 

Reference to Figure 1 shows that if the sample is left in an infinite bath such that 
_ the concentration just within the surface of the sample is constant at C,, the amount 
of penetrant in the sample reaches the value 4/C, when Dt//2=0-025. If the 
sample is insulated at this time the total content remains fixed and the surface 
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For the plane sheet excess absorption is 


the difference between /VW;¢/lc, at the 

time of transfer and at equilibrium. 
For the cylinder excess absorption is the 

difference between M;/7aC, at the 


plane sheet. 


The values of Dt/l? indicated are those at 
which the surface concentration falls 
to within 1% of its equilibrium value. 


time of transfer and at equilibrium. 
In each case M; is the content of the 

sample at time t. 
concentration falls according to equation (19). ‘The complete time variation of 
surface concentration is shown in Figure 6, for equilibrium contents 3/C,, 2/C,/3 
and 9/C,/10. The times taken for the surface concentration to fall to within 1% 
of its equilibrium value are indicated, and comparison with Figure 4 shows these 
times to be slightly greater than the corresponding times to condition to 99% 
equilibrium content by two-stage conditioning. Instead of insulating the sample 
the same result can be achieved by causing the concentration in the bath to vary 
in such a way that the surface concentration of the sample behaves acccording to 
the curves of Figure 6. Since the steady state is approached at approximately the 
sameé rate as in the two-stage conditioning, this method may in some cases prove a 
useful alternative to the two-stage method. 


§5. EFFECT. OF VARIABLE DIFFUSION OR THERMAL PROPERTIES 
The above mathematical solutions and discussion are based on the assumption 
that the diffusion coefficient, or thermal properties in the case of heat flow, are 
constant. There are many systems, however, for which this assumption is not 
justified, but where, for example, the diffusion coefficient varies considerably with 
concentration of penetrant. ‘This is frequently the case when water vapour cr 
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organic vapours or liquids diffuse into polymer sheets, though there is a lack of 
accurate measurements of the concentration dependence. When the diffusion 
coefficient is dependent on concentration it is no longer possible to write down 
formal mathematical solutions of the diffusion equation but a numerical solution 
can be obtained for each particular case. To illustrate in a general way the effect 
on the rate of conditioning of a diffusion coefficient which increases with increasing 
concentration of penetrant, a solution has been evaluated by a numerical process of 
integration (Crank and Nicolson 1947) for a diffusion coefficient D given by 
D=D,(1+10C/C,) where D, is constant and C and C have the same significance as 
in §4. In Figure 7 the curves for infinite bath and two-stage conditioning for this 
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variable D are compared with the corresponding curves for a constant coefficient 
D =D, taken from Figure 2, for the case in which the equilibrium content of the 
sample is 4/C,. Both methods lead to quicker conditioning for the variable 
diffusion coefficient than for the constant one, and furthermore, the advantage of 
using the two-stage method is slightly enhanced when D increases with increasing 
concentration of penetrant. Thus the ratio of the times of conditioning to the 99%, 
limit by the two methods is 5-14 when D= D,(1 4+ 10C/C,) and 5-00 when D= Dp. 
It can easily be seen from general reasoning that the effect of a diffusion coefficient 
increasing with concentration must be qualitatively in this direction, since the 
average value of D is greater than Dy and so both methods of conditioning must be 
quicker. Furthermore, the effect of using a high concentration in the first bath 
means that the relatively high values of D at high concentrations are allowed to 
influence the rate of conditioning, whereas in the infinite bath method concentra- 
tions above 4C, do not appear. In general therefore in a system in which the 
diffusion coefficient increases with increasing concentration the advantage gained 
by using two-stage conditioning is even greater than is indicated by the results for a 
constant coefficient presented in § 4. 

When D = D,(1 + 10C/C,) the excess absorption at the time of transfer does not 
differ greatly in the example considered from the corresponding excess when 
D=Dg, being about 0-092 in one case and 0-117 in the other. This suggests that 
in many cases, even though the diffusion coefficient is variable the mathematical 
solutions for a constant diffusion coefficient will be a useful guide to the amount 
of penetrant which the sample should contain when transferred from the first bath 
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) tothe second. In other cases, of course, it may be necessary to explore the system 
| to find the optimum time of transfer, either by carrying out the process of numerical 


¢ solution of the differential equation for the appropriate form of diffusion coefficient, 


| if this is known, or by practical trial and error. 


§6. CONDITIONING OF A CYLINDRICAL SAMPLE 


The four different methods of conditioning can of course be used if the sample 
approximates to the form of an infinitely long circular cylinder of radius a. The 
general reasoning is precisely as for the plane sheet and the appropriate formulae for 
the cylinder will be presented with a minimum of explanation. 

The basic equation for radial diffusion in a cylinder is 


ac _Da/, 2c 
Cf 7 Or Or 
In the infinite bath, taking Cy as before for the concentration of penetrant just 


within the surface of the cylinder, the required solution of (20) is (Carslaw and 
Jaeger 1947) 


= Jo(Ym7/@) | 
C=C {1-2 De exp(— De tle), — wases: 21 
e m=1 Gm? (Gm) Pt 2 Mf ( ) 


when the initial concentration throughout the cylinder is zero. q,,, is the mth root 
of the equation J,(q,,)=0. The integrated form of (21) is 


M,=7a°C, {1- as ae (- Dai’). Ree (22) 
m=1%m 


The first part of two-stage conditioning is governed by equations (21), (22) with C, 
substituted for Cy if C, is the concentration of penetrant at the surface of the 
cylinder in the first bath. Diffusion in the second bath is subject to the conditions 


C=C, =constant, T=, fT ae We as oa (23) 

‘C= 5/Cay O'=cr <a; Pe ee Ve (24) 
where ¢, is the time of transfer from the first bath to the second, and in general the: 
solution of (20) with (23), (24) is 
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where J,(g,,) =0 now. 
For the particular case where f(r) is given by substituting tf =4,, Cy =C;, in (21), 
we find 


[fer Volar’ Ia) dr’ = Cal —exp(— Digital®) W(Gu) Ida (26) 
and (25) becomes 
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of which the integrated form is 


M,=7aC, {1-4 3 Zexp{—Dai(t -1)/4°}} 


n=1%4n 
+4nCya? & = exp{ —Dq2(t—t,)/a?}{1 —exp(—Dq>t/a")}. .....- (27 a) 
n=14n 


For the finite bath problem we use cylindrical polar coordinates and assume the 
cylinder to occupy the space r<a, while the cross-section of the bath in which it is 
immersed is A, excluding the space occupied by the cylinder. ‘The concentration 
in the cylinder is initially zero and in the bath is always uniform. ‘The equation 


expressing the conservation of penetrant is 
= +2n | rCdr= — gr 0 ee eee (28) 


where C, is the initial concentration just within the cylinder surface and is the 
partition factor for the penetrant in the sample and the bath as before. he 
required solution of (20) for these conditions is 


4(1 +a) exp (= Dgnt/a*) Jo, at) 
C= Goa x a Rh ee ie hak 9A 8d | 29 
02 ee 
where the q,,’s are the positive roots other than zero of «q,,J9(q,,) +2/4(q,,) =0, and 
C ,,=aC,/(1+«), «=A/(za*k). The integrated form of (29) is (Wilson 1948) 
M, S 4a(1 +a) 
Ne na14+4a+a2q2 
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where M,,=C,A/R(1+«). ‘The solution of (20) for a cylinder insulated at time t, 
and with the initial condition C=f(r), O<r<a, t=t,, is (Carslaw and Jaeger 
1947) 


C= a]. r' f(r’) dr’ i z coy — Dp*(t —ty)/a?} ge fs flee bar! la) ar'} , 


where the p,,s are the roots of 


SC ees a Vt (32) 


If, prior to insulation, the concentration just within the surface of the cylinder is C), 
so that f(r) is given by evaluation of (21) with C)=C, and t=#,, then it follows 
readily that 


| rf") dr = 3Cya2}1—4 x Ce Daitsla®)) Mg: (33) 
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using the relation if rr 1J(Br) dr =p" +17, _,(Br). Also by using the relations 


and the fact that, with (32), | rJ ((p,1/a) dr =0, it can easily be shown that 
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so that finally the concentration at the surface of the cylinder after insulation is 
given by 


Be igh 
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By numerical Grae of the mathematical solutions of this section the four 
methods of conditioning a cylindrical sample may be compared. The general 
behaviour of the overall rate of absorption and concentration—distance curves is 
closely similar to that of the corresponding curves for the plane sheet. For any 
given equilibrium content of the cylinder there is an optimum point of transfer 
from the first bath to the second in two-stage conditioning, which leads to a shorter 
conditioning time than that needed for finite bath conditioning, this latter being 
quicker than conditioning in an infinite bath. In insulated conditioning the 
surface concentration falls to within 1°% of its equilibrium value in a time which is a 
little greater than the time of conditioning to 99° equilibrium content in two 
stages. The main features of practical interest for the cylindical sample are 
summarized in Figures 5 and 8. Figure 8 shows how the time of transfer and the 
time to condition to 99°% of the equilibrium content by two-stage conditioning 
depend on the equilibrium content required, the concentration in the first bath 
being fixed. ‘The upper curve in Figure 5 shows the excess absorption at the time 
of transfer as a function of equilibrium content. For asystem in which the diffu- 
sion coefficient or thermal properties are constant the data of Figures 5 and 8 indicate 
how to carry out two-stage conditioning of a cylinder in order to achieve the 
optimum conditioning time. 


Wee IVa LON S ON HE PRACTICAL APPEICATTON OF 
THE MATHEMATICAL TREATMENT 

All the mathematical solutions developed in this paper refer to the conditioning 
of a sample which is initially free from penetrant. The treatment could clearly be 
generalized, however, to cover cases in which the initial content is not zero. 

A more serious restriction on the type of system to which the theoretical pre- 
dictions may reasonably be expected to apply, arises out of the implicit assumption 
that there is no hysteresis in the vapour Pressure —content isotherm. Thus 
comparing for example “infinite bath conditioning” with “two-stage condition- 
ing’’, the former process is one of absorption only, whereas the latter is one of 
absorption followed by desorption. If the material shows hysteresis, the actual 
contents at the end of the two processes will be different, even though the vapour 
pressures are the same. These considerations limit the direct application of the 
treatment to systems in which hysteresis is not an important factor. 

Somewhat similar restrictions will apply to the two-stage heat conditioning 
process if the material undergoes, at high temperatures, any physical change 
which is not entirely removed as the temperature falls again. 


REFERENCES 
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A Travelling Wave Linear Accelerator with R.F. Power 
Feedback, and an Observation of R.F. Absorption by Gas 
in Presence of a Magnetic Field 


A travelling wave linear accelerator usually consists of a high frequency power source: 
feeding power into a waveguide accelerating tube which is terminated by a dummy load’ 
to absorb without reflection the power reaching the output end. It has been shown 
(Harvie 1948) that there is an optimum length of accelerating tube, which corresponds 
to an attenuation of approximately 1:25 nepers, but unfortunately there are sometimes 
practical difficulties in making accelerators as long as the optimum length. This note 
describes means for overcoming some of these difficulties, in which an accelerator of less. 
than optimum length retains its efficiency by utilizing the power which is usually wasted 
in the dummy load. 

The system is shown diagrammatically in Figure 1, in which arrows indicate the 
direction of power flow. Power is fed back from the output end of the accelerator and. 
combined with the power from the source in a suitable bridge circuit. It is sometimes 


Thermocouple Power Monitor 
[) ypVariable Probe Magnetron 


Phase Shifter 
Steel Load 


Rat-race 
Brid Feedback Line (n=l) 
nieee Phase Shifter (adjustable in 
operation): 
Pre-set 
7 Phase Shifter 
HF Power Source (a) Travelling Wave 
W,=0 Accelerating Waveguide AQutput 
‘ End 
Absorbing Load 4mev Linear Accelerator 
Figure 1. Diagrammatic arrangement of Figure 2. Feedback system applied to 4 Mev. 
linear accelerator with feedback. linear accelerator. 


desirable that arms (1) and (2) of the bridge should be mutually conjugate, in order to make 
the impedance presented to the source constant during a transient building-up period and 
independent of the accelerator characteristics. It can then be shown that the bridge 
circuit must include a dissipative element which is shown connected. to arm (3) externally; 
under steady state conditions no power is dissipated in this element. Steady state power 
fluxes are shown on the diagram; the power flux entering the accelerator (1+-n”)Wg, 
is always greater than the power supplied by the source and 7 is determined by the attenua-- 
tion of the accelerator. Calculations show that the “ efficiency’’ of the accelerator, 
V?/WL, increases as the length is diminished and is slightly greater than that of an optimum: 
length accelerator with a dummy load at the end. 

Bridges of the type required have been studied previously (Tyrell 1947, Smullin and 
Montgomery 1948) for the special case where n=1, and examples will be familiar under 
the names of “ magic T’” and “ Rat-race’’. There seems to be no great difficulty in 
designing bridges along similar lines for other values of n. 

In a pulse operated system it is necessary for waves to make several transits of the 
accelerator and feedback loop, after application of the pulse, before steady state conditions 
are approached. Calculations show that the build-up time is independent of the length 
and is approximately the same as in a conventional accelerator of optimum length. | If’ 
the frequency of the source is varied the power fed back from the end does not return in 
the correct phase and the circulating power is reduced. It can be shown that the permissible 
frequency bandwidth is again approximately the same as in an ordinary accelerator of’ 
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optimum length; however, by adding an adjustable phase changer in the feedback path 
it is possible to maintain the correct circulating power at any frequency and the bandwidth 
is then only limited by the wave-to-electron phasing in the usual way. 

An advantage of the system is that the accelerator can be made appreciably shorter~ 
than the optimum length without losing efficiency, which leads to increased dimensional 
tolerances on the waveguide. A shorter waveguide also has a wider frequency bandwidth, 
as far as phasing of the wave to electrons is concerned, and it may also be an advantage 
in practice to allow the frequency of the source to vary as it will within this band while 
the operating conditions are kept right by adjusting the phase changer in the feedback 
loop. 

For the case of a bridge circuit designed for n=1, theoretical curves can be plotted 
for the power in the various arms of the bridge when the phase and attenuation in the 
accelerator loop are varied. Experimental confirmation has been obtained at low c.w. 
power with a mock-up system on a wavelength of 3-2 cm., using a normal n==1 “‘ Rat-race ”’ 
and about 30 wavelengths of rectangular waveguide to represent in some measure the 
corrugated waveguide. A method has also been worked out theoretically and confirmed 
experimentally for obtaining a suitable signal which will reset the phase shifter in the 
feedback arm when the magnetron frequency changes. 

It has also been possible to try out the n=1 system on the 4 Mev. linear accelerator 
(Fry et al. 1948) since the r.F. loss in the accelerator is about 2:5 db. and can be made up 
to 3 db. with beam loading. The experiment had unfortunately to be carried out in a very 
short time and the external waveguide components were not matched as well as they should 
have been; consequently the thermocouple power monitors could not be relied upon to 
within 20°, although they were partly compensated for by the presence of standing waves. . 
It was not therefore possible to verify the power fluxes in all parts of the system. However, 
independent measurements showed that the accelerator operated as if it had an input 
R.F. power flux of 2:1-2:2 mw. when the magnetron power was only about 1-4mw. ‘The 
energy spectrum was identical to that previously obtained with the same input power 
but without feedback. The maximum power which could be obtained in the accelerator 
was 2:4 mw. and, with a somewhat lower frequency than normal, the mean energy was 
raised to 4:5 Mev. but with more than twice the normal spectrum width. ‘These experiments 
indicated that a well matched system would work as predicted, but with one reservation— 
that a somewhat better vacuum is required than for an equivalent accelerator without 
feedback. 

It was noticed that the loss in the corrugated waveguide soon after applying the R.F. 
power but without the gun in operation was very greatly increased when the focusing field 
was switched on, and moreover the effect happened quite critically at about one-third of- 
the normal focusing field and was not confined to any particular spot along the length of 
the corrugated waveguide. It has been checked that the power leaving the accelerator 
falls in the same manner for the system without feedback, but since there is no reaction on 
the input power its effect on this particular accelerator is not serious and had therefore 
not been observed before. The effect may be due to the presence of positive ions or electrons 
which spiral under the influence of the focusing field and consequently abstract very 
much greater quantities of energy from the R.F. wave than if they drifted straight to the 
walls of the tube. As the vacuum cleans up the effect entirely disappears and the operation 
is as predicted. 

This work was carried out as part of the programme of the Accelerator Group of the 
Atomic Energy Research Establishment at the Telecommunications Research Establishment, 
Malvern, and acknowledgment is made to the Director, A.E.R.E., for permission to: 
publish the results. 
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The Dielectric Properties of Ice at 1:25 cm. Wavelength 


It is known that ice exhibits a loss-maximum dependent upon temperature, in the 
frequency range 104 to 10° c/s. (Granier 1924, Errera 1924, Smyth and Hitchcock 1932, 
Eder 1947): it has also been shown by Lamb (1946 a) that, at —5° c., the resulting absorp- 
tion is characterized by a single relaxation time. Lamb (1946a) further measured the 
‘dielectric properties at a free space wavelength of 3 cm. over the temperature range 0 to 
—50° c., using an Ho, cylindrical cavity type of resonator. Under these conditions the 
permittivity was sensibly constant and was given as 3:0;; further experience with the 
-apparatus used in these experiments, however, has indicated that this value should be 
increased by approximately 4°, owing to the use of a metal cup to contain the ice specimen 
within the resonator. The corrected permittivity value of 3-1, is appreciably greater 
than the square of the optical refractive index, x?=1-72 (International Critical Tables I, 
166), and the measurements have now been extended to a free space wavelength, A=1:25 cm., 
and to temperatures as low as —195°c., in an attempt to observe a second dispersion. 
An Ho, cylindrical cavity resonator has also been used in these experiments (Lamb 1946 b), 
but in this case the ice specimen was formed to fit loosely in the cavity and rested on the 
movable piston. Each specimen was prepared from water which had been distilled once. 
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Curve 1. Measured permittivity values at A=1°25 cm., 
without correction for change in specimen 
thickness with temperature. 

Curve 2. Measured tan 6 values at A=1°25 cm 

Curve 3. Tan 6-temperature curve for A=3 cm. repro- 
duced from Lamb’s paper. 


The results obtained are given in the Figure, in which the tan 6-temperature curve 

or A==3 cm. is reproduced. At any given temperature the loss-factor value at A=1:25 cm. 

‘is slightly less than that at the higher wavelength, but its variation with temperature is of 
a similar form. No appreciable change in permittivity was observed over the temperature 
range covered by the experiments : the small decrease at lower temperatures shown in the 
Figure could be accounted for by the change in specimen thickness with temperature 


for which no correction was applied. The permittivity near 0°c. is 3-18 (12%), in 
agreement with the corrected value for A=3 cm. 
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A decrease in temperature would displace a dispersion due to dipole relaxation towards 
a region of longer wavelength. Since, however, there is no evidence of a decrease in 
permittivity (nor of an increase in tan 5) at A=1:25 cm., even for temperatures as low as. 
—195°c., it is unlikely that the difference between the measured permittivity and the 
square of the optical refractive index can be explained by a further dipole relaxation. 
The difference could be attributed to a dispersion caused by vibration of the molecule 
about its equilibrium position in the crystal lattice at a resonance frequency considerably 
greater than 24,000 Mc/s. (A=1-:25 cm.): an alternative explanation would be a com- 
paratively large “ atomic’”’ polarization. ‘The contribution to permittivity at A=1-:25 cm. 
arising from either of these latter effects would be sensibly independent of temperature. 

One of the authors (A. T.) wishes to thank the British Electrical and Allied Industries 
Research Association by whom he was employed during the progress of this work. 
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REVIEWS OF BOOKS 


Dipole Moments, by R. J. W. Le Fiver. Pp. 117. Second Edition, revised and 
reset. (London: Methuen, 1948.) 5s. 


Since the first appearance, ten years ago, of this excellent little book, the reviewer has 
consistently recommended it to and for honours students in chemistry. Seeing that the 
second edition differs hardly at all from the first, he will continue to do so. The only. 
changes worthy of note are the following. 

At the very beginning of the*ook there is a new short section on the relationship of 
dielectric constant to other properties. This contains a sketchy account of certain empirical 
relations, so ancient as to be of scant interest. The reviewer cannot help wondering 
whether this insertion was by request. 

The very useful chapter on experimental methods of measurement has been made yet 
more useful by the inclusion of further practical information. 

The chapter on molecular structure contains new data on several pairs of geometrical 
isomers which have been investigated by the author and collaborators. "The author makes 
the interesting suggestion that the rate of inter-conversion of geometrical isomers in 
solution could be determined by measuring the dielectric constant. 

The reviewer’s only criticisms are directed not at this particular book, but against 
certain features common to the whole literature of this subject. 

One does not measure a dipole moment. One measures dielectric constants and 
refractive indices. "Then by unnecessarily cumbrous algebraic formulae one calculates a 
dipole moment. 

If the dipole moment thus calculated is found to differ according to the solvent, the 
cause may well be neither solute—solute interaction nor solute—solvent interaction, but 
just the inaccuracy of Debye’s approximations, whose main justification is the lack of any 
better. This inadequacy was realized and stated by Debye, but seems to be forgotten or: 
ignored by the majority of those who use his formulae. 

When Onsager developed an improved theory of the dielectric constant of pure liquids,. 
he indicated briefly that the treatment might be extended to solutions, but no one has yet 
followed this lead. We still lack a satisfactory theory of the dielectric constant of liquid. 
mixtures. 1s WES. 


274 Reviews of Books 


Progress in the Theory of the Physical Properties of Glass, by J. M. STEVELS. 
Pp. xi+104. (New York, Amsterdam, London, Brussels: Elsevier “4 
Publishing Company, Inc., 1948.) 10s. net. y 


Theories of the physical properties of glasses were almost non-existent until the pioneer 
work of Zachariasen and of Warren during the eight years immediately before the war. 
The present book describes the advances made in Holland during the war mainly through _ 
the efforts of Stevels and his associates at the Philips Research Laboratory, Eindhoven. — 
In his preface the author makes it clear that no attempt has been made to take into account 
all the American and British literature published since 1939, although some additions from — 

‘these sources were made during the correction of the proofs. 

The first chapter presents a general account of glass in the solid state based mainly on — 
the experimental work of Warren and the ideas of Zachariasen, who stated rules for deciding 
which inorganic oxides should form glasses in cases where the structures of the corresponding 

crystalline forms are known. The author’s claim that complete agreement with the 

Zachariasen theory is found is of doubtful accuracy; the Zachariasen theory should be 
regarded rather as an excellent guide to probable glass-forming substances, but not yet 
-completely proved or entirely free from ambiguity. 

A long chapter on the density of glass is based on the theory of Stevels that the volume 
of a block of glass containing one gram atom of oxygen ions is in many cases given by an 
empirical formula involving the number of oxygen ions divided by the number of network- 

»forming ions (Si*+, B** etc.) and a constant depending upon the nature of the network 
formers, but to a first approximation not upon the nature and the amount of the network 
modifiers (Nat, Ca?+ etc.). The significance of the formula in relation to structure is 
discussed at length for silicate, borate, borosilicate, boroaluminate, phosphate and germanate 
glasses. The chapter forms an impressive example of the value of a semi-empirical attack 

-on the physical properties in advancing our knowledge of glass structure, even though 
neither the formula nor the values of the constants can be derived theoretically. 

Chapter 3 deals in a geneial way with the electrical conductivity of glasses and provides 
a derivation of the Rasch and Hinrichsen law. Chapter 4 discusses dielectric losses in glasses ; 
the theoretical background provided is the recent one of Gevers and du Pré, together with a 
theory of the relation between loss angle and glass composition. The theories should 
stimulate further research, particularly on the losses at low temperatures and over greater 
‘ranges of frequency or glass composition. 

The fifth chapter deals with molecular refraction, the refraction being regarded as the 
sum of the contribution of each ionic constituent. The’theory is restricted to certain ranges ~ 
of composition only and, although the ionic refractivity of silicon is taken from the previous 
work of Fajans and Joos, the values for the oxygen, sodium, calcium and potassium ions 
are adapted to fit the experimental results. 

It is easy to be critical of newly developed theories, and in the case of glass particularly 
so, in that so little sound theoretical work has been done in spite of a great number of 
technological measurements of physical properties. Dr. Stevels is to be congratulated 
therefore on his efforts to develop the theoretical aspects of glass physics. "The book 

-deserves careful study by all interested in glass science, and will serve its purpose if the 
limitations of the theories he puts forward stimulate further investigation. 

J. E. STANWORTH,. 
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The Performance of Infra-Red Spectrometers as Limited by Detector Characteristics, 
by E. F. Daty and G. B. B. M. SuTHERLAND. 


ABSTRACT. The performance of the average infra-red spectrometer is conditioned 
mainly by the characteristics of the detector employed. Assuming a minimum usable 
energy level at the exit slit, the square of the minimum frequency band embraced by the 
exit slit is proportional to the minimum energy level and to the focal length of the collimator, 
and inversely proportional to the angular dispersion of the prism or grating, the area of 
the collimator mirror, the energy flux at the entrance slit and the slit length. If, however, 
the time of response of the detector is taken into account, the inter-relation of spectrum 
scanning speed, signal-to-noise ratio and resolving power must be considered. For this 
purpose a simplified model of an infra-red detector has been used, defined by three parameters. 
This allows the relative merits of detectors to be assessed and the performance of any 
‘given spectrometer with a specified detector to be estimated. ‘The inter-relation of resolving 
power and spectral scanning speed has been examined for “‘ chopped-beam ”’ and “ video ”’ 
presentation. - 


The Luminescence of Wetted Solids, by J. Ewes and G. C. FARNELL. 


ABSTRACT. An earlier observation that the association of a non-luminescent pure 
liquid with a non-luminescent solid produces fluorescence and phosphorescence in the 
liquid has been confirmed. ‘The luminescent spectra of water, heavy water, and alcohols, 
as well as some other liquids with different negative radicals when adsorbed by various 
finely divided ionic solids have been examined and recorded. ‘The positions of the maxima 
of the main broad band are very similar for all the hydroxy liquids even with different solids 
(band peaks in the range 4300-4500 a.) and are different for the same solid with liquids 
having different negative radicals. Preliminary experiments on the relation between bright- 
ness and relative amount of liquid indicate that the phenomena are associated with the 
active sites on the solid surface, and afford evidence of the existence of these and of their 
.differences on the same solid. 

Excitation spectra confirm preliminary observations that the adsorbed molecules have 
very abnormal absorption regions indicating considerable changes in energy levels. 'T'wo 
‘suggestions are made to explain the phenomena: (a) the energy changes may be correlated 
with the heat of adsorption on the active sites; (b) the adsorbed molecules are bonded 
to the solid by the positive radicals leaving an almost free negative ion. 
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The Resolving Power and Intensity Relationships of the Fabry Perot Interfero-— 


meter with Silvered Reflecting Surfaces, by R. J. Bricut, D. A. JACKSON 
and H. Kuun. 


ABSTRACT. By the use of a Fabry Perot etalon with very small spacing, the widths of the 

interference fringes were measured under conditions where the light could be considered as. 
practically monochromatic. In this way, the resolving power which can be achieved with 
silver deposits of known optical density was determined for light of five different wavelengths. 
between 6438 and 36104. The results can be used for selecting the most favourable 
thickness of silvering of etalon plates, namely that which gives the resolving power required 
for a given purpose without unnecessary loss of intensity. The application of the results to 
multiple etalons is discussed. 


Intensity Distribution for Bands of the y-System of MgH*, by M. E. PILLow. 


ABSTRACT. ‘The intensity distribution for bands of the y-system of MgH* is deduced. 
theoretically from the wave functions for the various energy states, and is compared with the 
distribution of the observed bands. 

Lines of the (2,3) band, not previously noted, have been observed. 


Nuclear Disintegrations in Photographic Plates exposed to Cosmic Rays under Lead — 


Absorbers, by E. P. GeorceE and A. C. Jason. 


ABSTRACT. Ilford Nuclear Research photographic plates have been exposed at the 
Jungfraujoch, altitude 3,457 metres, under lead screens of various thicknesses up to 28 cm. 
The intensity of the stars, and its variation with the number of tracks forming them, have 
been recorded for each thickness of absorber. ‘The absorption curve of the star intensity 
shows that the range in lead of the radiation causing the stars is 310-+-20 gm/cm?. ‘The size 
distribution is apparently the same under all absorbers. 

The intensity and size distribution have also been recorded in similar plates exposed at 
sea level. ‘The size distribution is again the same as that found for the sbigh altitude plates, 
but the range in air of the primaries is 150-+-7 gm/cm?. 

These results are consistent with the generally accepted view that the majority of the 
stars is caused by neutrons. 


Experiments with the Delayed Coincidence Method, including a Search for Short- 
lived Nuclear Isomers, by D. E. Bunyan, A. LuNpBy and D. WaLKER. 


ABSTRACT. "he delayed coincidence method has been studied in some detail and. 
employed in a search for short-lived isomers in the product nuclei of B-emitters. The 
influence of counter time-lag fluctuations, which determine the short half-life limit in this. 
method, is discussed in relation to the “‘ differential? procedure for recording the delayed 
coincidences. It is shown that the time lags exhibited by a Geiger-Miiller counter are 
distributed approximately according to a Gaussian law. A differential delayed coincidence 
recorder has been constructed in which the overall time-lags in the two input channels due 
to “‘ rise-time ”’ effects have been stabilized. At the coincidence stage of this recorder the 
quadratic mean deviation in the time-lag difference between the two channels due to 
“rise-time ”’ effects has been kept as low as ~5X10-*®sec. Details are given of the 
experimental procedure used in searching for short-lived isomers. The range of possible 
half-lives most thoroughly investigated extends from about 3 x 1077 sec. to about 10-3 sec. 
for isomeric transitions giving conversion electron energies greater than about 100 kev. 
Possible improvements in the range of the experiments are discussed. The short-lived 
isomers of 1*”T'a and 18*Re have been confirmed. The half-life of the metastable state in 
187Re has been found to be (5:26+0-12) x 10-7 sec. Negative results are reported for 19 
isotopes. In particular the metastable states reported by other authors in !%Hg, !41Pr, and 
124Te are not confirmed. 
The half-life of ThC’ has been found to be (3:04-+0-04) x 10~? sec. 
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